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Introduction 


y  This  report  deals  with  the  mathematical  aspects  required  for  the 
establishment  of  a  computer  program  able  to  calculate  all  stresses 
and  displacements  in  a  four  layered  system. 

The  materials  of  the  different  layers  may  be  isotropic  or  cross - 
anisotropic. 

The  interface  conditions  cover  all  the  cases  from  full  friction 
to  full  slip  included  partial  friction. 

The  bottom  of  the  fourth  layer  is  considered  to  be  fixed  (no  vertical 
deflections).  'i 

( / 

^The  report  is  based  on  A  r  f:  w  ( v  _ ^ 

^  existing  material  :^"otropic  multi layer^wr)^*^RPISTEfrr^943)  ~ 
^ ^  ^and~an~isotropic  muTtilayer  theory; ^VAW--^UWEtAEftT-,-^:9^  < 

fy  original  research  interface  conditions  (fixed  bottom, 

partial  friction)  and  satisfactory  convergency,  thus  complete 
accuracy,  at  the  surface  and  in  the  first  layer  of  the  system. 


^  r-f  V /pus 


s-^a  rcji 


This  report  contains  ^htee  partS^  C^) 

Rart  rn  a  theoretical  outline  whe^n  the  basic  equations  are  given, 
the  specific  boundary  conditions  discussed  and  the  particu¬ 
lar  numerical  problems  relied  to  the  accuracy  at  the  sur¬ 
face  and  in  the  first  layer>v solved.  (  2.) 

/  r 


Part  Z 

C 


the  general  mathematical  analysis  of  the  chosen  numerical 
solution  and  a  description  of  the  programs  and  thsir 
utilization.  O.  ^  ^  ^ 

S  - 


Pa/t  3  : ^(appendices)  the  detailed  mathematical  and  algebraical  analy^f-5 
^-frisses  for  the.  different  considered  cases, isotropic  or 
anisotropic,  full  slip  or  not.  _ 

The  programs  are  written  in  FORTRAN  77  and  run  on  IBM  PC  or  (aM^her 
compatible  equipment,  /fr’v  ;  Voty^n^'r.  4  cJe  Jrf  ^  — 
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The  basic  equations  and  fundamental  hypothesisses. 
The  boundary  conditions 

1.2.1.  The  partial  friction  condition 

1.2.2.  The  fixed  bottom  condition 

Solution  of  particular  numerical  problems 

1.3.1.  The  full  slip  condition 

1.3.2.  Over  and  under  flow  problems 

1.3.3.  The  vertical  deflection  at  the  surface 

1.3.4.  Accuracy  in  the  first  layer 


PART  2  -  Numerical  resolution 

2.1.  The  mathematical  analysis 

2.1.1.  Full  or  partial  friction  at  the  interfaces 

2.1.2.  Full  slip  at  the  first  two  interfaces 

2.1.3.  The  detailed  mathematical  analysis 

2.2.  The  numerical  procedure 

2.3.  Flow  sheet 

2.4.  Particular  features 

2.4.1.  Verification  of  the  accuracy  of  the  results 

2.4.2.  Converoency  of  the  Bessel  functions 

2.4.3.  Underflow  problems 

2.5.  Use  of  the  programs 

2.5.1.  Generalities 

2.5.2.  Source  version 

2.5.3.  Running  of  the  programs 

2.5.4.  Miscellaneous 

2.5.5.  Values  for  the  degree  of  anisotropy  and  the  partial  friction 
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APPENDIX 


Analysis  of  the  isotropic  case  with  full  and  partial  friction 

Analysis  of  the  isotropic  case  with  full  slip 

Analysis  of  the  anisotropic  case  with  full  and  partial  friction 

Analysis  of  the  anisotropic  case  with  full  slip 

Explanatory  notice. 


-  C^m  (2  -  4y^  -  mzle”’^  -  V^m  (2  -  4vi^  +  mz)e  ]  dm 


1+Mt  f  Ji  (mA)  .Ji  (nia)  2„mz  ,  „ 

U  --  -  J  ^  *  ^-i 


2  -mz 
me 


+  Cj-m  (1  +  mz)e™^  -  f)^m  (1  -  mz)e  1  dm 


where 


s  the  radius  of  a  circular  uniformly  distributed  load 
s  the  value  of  the  vertical  pressure 

s  the  horizontal  distance  from  the  axis  in  a  cylindrical  coordinate 
s  the  depth 
s  the  vertical  stress 

s  the  horizontal  radial  stress,  o«  is  the  circumferential  stress 


is  the  circumferential  stress 


s  the  shearstress 
s  the  vertical  deflection 


a  is  the  radial  (horizontal)  displacement 

El  is  the  Young  modulus  of  the  concerned  layer 
is  the  Poisson's  ratio  of  the  concerned  layer 
A^,P(  are  unknown  parameters  to  be  determited  by  the  boundary  conditions 
Jo  is  the  Bessel  function  of  the  first  kind  of  order  zero 

Ji  is  the  Bessel  function  of  the  first  kind  of  order  one 

m  is  an  integrating  parameter 


In  the  case  of  a  cross  anisotropic  body,  they  are  obtained  from 

(VAN  CAUWFLAERT,  1983)  : 

m 

^  .  pa  I  rim 


This  stress  function  differs  fundamentally  from  the  proceeding  one  : 
indeed  in  putting  .6^  =  1  in  it,  we  do  not  obtain  the  stress  function  for 
the  isotropic  case.  We  conclude  that  the  tv^o  cases  must  be  handled  in  a 
separate  way. 


The  stresses  and  displacements  are  given  by 
O2  =  pa  I  Jo  (nw) . Ji  (ma)  [n^d  + 

O 

+  +  p^l  ]  dm 


pci  Jo  (mA) .  Ji  (ma)  (n^ll  +  +  8^m*e  ) 

^  O 

k./Ih;  -  U^jI  i  2  ^;fnz  ,  f,  .  2  ,  , 

t —1 — i - '—L.  (C-i-ni  e  ^  +  P:i  m  e  -c  )  ]  dm 

-  Pt 


.  f  Ji  (m-l )  .  Ji  (m^l  (1  .  I  f  1 

J  -- — +  I't)  lAi 


2  niz  .  D  2  -mz 
m  f  +  8^-  m  c 


+  +  V^i^ni^e  -c*”  )  dm 


pa  [  Jo  (frw) .  Ji  (ma)  (  C;  +  P.- 4- 

m 

pa  j  JilwUli.lrol  ,  *  B,„>e-"’Z  ,  0. 


+  P^  6^m^e  1  ( 1+p^  I  dm 


^Az  =  -  P«  j  Ji  I'M)  .J,  (ma)  ln^(l  +  y^l  (A^m^e'"^  -  B^m^c“'"^) 


+  pj  (C^  -  P^  S^m*e''^-i'"^n  dm 


.  Liji'  pa  f  M»«!.^!!4l  „^.(l  t  „,| 

'  o 

,  (”<  +P<)  ,  7.  i;mz  „  ,  -j.  mz,,  , 

+  - - ^ -  (f/<4;m  c<  -  P;i:m  e  <  1  dm 


I  'im 


•'.  -\V, . 


'«'*  %'  V  ’ 

'./.■•  w 


6. 


u  = 


(1  +  Rt)  +  M/) 


J.lwl.J,li».|  I  ■  ^  B, 

J  m  ^ 


2  -mz 
;m  e 


,  r>  2  i;  mz  .  n  ,  2  -  i;  mz  ,  , 

4  r^-A^m  (>  <  +  ('  ^  1  nm 


where 

is  the  degree  of  anisotropy,  the  ratio  between  the  vertical 
and  the  horizontal  Young  moduli  of  the  concerned  layer  ; 

is  Poisson's  ratio  expressing  a  strain  in  the  horizontal 
plane  induced  by  a  stress  in  the  vertical  direction  ; 

is  the  index  of  anisotropy. 

The  anisotropic  relations  arc  established  in  tlie  assumption  tliat 

-  the  shear  modulus  in  the  vertical  plane,  is  related 

to  the  other  elastic  constants  by  (BAPf)rN,1963  ; 

VAN  (:AUWELAFRT,19P3)  : 

1  -  1  +  •»  ^P<' 

GKzi  El 

-  Poisson's  ratio  in  the  horizontal  plane,  v,  is  related 

to  Poisson's  ratio  in  the  vertical  plane  by 

(rFTlMIF,1973)  : 

V,-  = 


thl 


ht 


W  I 


THE  BOUNDARY  CONOITIONS 


1.2.1.  The  partial  friction  condition 

Let  us  consider  an  layered  system,  consisting  in  (n-1)  layers  of  a  finite 
thickness  built  on  a  semi-infinite  body. 

For  each  layer  exists  a  stress  function  ‘Jk' ( A/ Q P,- )  v/ith  4  unknown  para¬ 
meters  :  the  total  of  unknown  parameters  is  4n. 

Two  parameters  depend  on  the  shape  of  the  load  at  the  surface 

oz  =  /  (p)  for  z  a 

lAz  -  0 

At  infinite  depth  stresses  and  displacements  must  vanish  and  thus  A„  and 

=  0. 

We  remain  with  4  (n-1)  parameters  to  be  determined  with  4  conditions 
at  each  interface. 

The  hypothesis  is  introduced  at  this  stage  that  under  effect  of  tlie  load, 
the  layers  remain  individually  fully  in  contact,  what  we  express  by  imponing 
that  at  the  bottom  of  each  layer  and  at  the  surface  of  next  layer  vertical 
stresses  (oj),  shearstresses  i'^nz)  vertical  displacements  (u')  are  identic 
The  fourth  interface  condition  depends  on  the  relative  adhesion  in  the 
horizontal  plane  between  the  considered  layers. 

The  two  extremes  are 

-  full  continuity,  expressed  by  setting  that  the  horizontal  displace¬ 
ment  (u)  are  identic  ; 

-  frictionless  interface  by  considering  the  interface  as  a  principal 
plane  and  thus  by  setting  the  shearstresses  equal  zero. 

Partial  adhesion  has  been  temptatively  introduced  by  several  authors, 
utilizing,  in  the  same  v/ay  as  WESTERGAARD  (1926),  a  relation  betv/een 
horizontal  displacements  and  shearstress  : 

K  (ii(  -  U{  +  j}  =  Taz^' 

where  u(  is  the  horizontal  displacement  at  the  bottom  of  the  i-th  layer 
and  (i(+j  that  at  the  surface  of  the  (i+l)-th  layer. 

We  shall  prove  that  such  a  relation  cannot  be  correct  in  the  case  of  a 
mul ti layer. 


One  has,  for  an  isotropic  body,  followinn  relations  betv/een  displacements 
shearstrains  and  shearstresses  : 


^  hi  =  y^:zl  .-  (  2(1+  p^  )/F^l 


3z 


<lz  ^ 


+  j  =  I  2  (1  +  ii,'+i  )/f  <7(  1  T^Zj  +  i 


We  know  from  the  interface  boundary  conditions  that 


W,-  =  10 


^fiZ(  ~  ^hzi+\ 


and  thus  that 


3  I  l  o  r  l+bi.  1+P/  +  1  • 

■  “«  +  i'  "  ^  ^  ■  Fy+1  ' 


l+Pc  l-*lJ<  +  i 


I  \  O  f  11  C  T  I  , 

-  “c+i'  =  ?  f  T~- - r —  1 


T/izv  •  ^z 


On  the  other  hand 
K  (U{  -  U(  +  i,  =  ’’^Az^ 


so  that 

T/tz^  •  dz  = 

which  solution  is 


2  K  I 


UUt 

E; 


UU.(>  1  ,  ‘ 

r  i 


=  exp  (  2K.z( 


l+V't 


Ubct  1 


)  1  .  /  (A) 


Comparing  this  solution  with  the  relation  above  for  the  shenrstress, 
v/e  conclude  that  the  obtained  expression  cannot  be  deduced  from  a  stress- 
function  that  is  still  solution  of  the  compatibility  equations. 
Compatibility  is  thus  nos  respected  and  relation  K{u(  - 
cannot  be  accepted. 


9. 


Our  meaning  is  that  the  only  v/ay  to  express  partial  adhesion  consists  in 
writing 

H  *  <“4+  1 

with  K  e  1  o I 

When  K  ^  I,  one  has  full  continuity 

1^/1,  one  has  partial  continuity  (normally  K>  1) 

Full  slip  is  then  a  completely  other  case,  for  which  another  program  has 
to  be  written,  what  the  before  mcntionned  authors  tried  to  avoid  by  intro¬ 
ducing  their  particular  relation. 


1.2.2.  The  fixed  bottom  condition 

The  boundary  conditions  discussed  in  previous  alinea  implicate  that  the 
last  layer  of  the  multilayer  is  considered  as  a  semi-infinite  body. 

One  can  also  consider  the  case  of  a  multilayer  built  on  an  uiicic'fomiable 
body,  that  thus  any  vertical  displacement  vanishes  at  the  contact  face 
with  the  undeformable  body  :  we  shall  call  this  a  fixed  bottom  condition. 
This  condition  can  be  introduced  in  several  manners  and  thus  demand  a 
detailed  analysis. 

A  vertical  displacement  is  obtained  by  integration  of  the  vertical  strain  : 
w  •  I  Cz  .  dj 

It  would  not  be  correct  to  resort  to  an  integration  between  limits, 
such  as 


where  H  could,  for  example,  be  the  depth  at  which  we  want  the  bottom  to  be 
fixed.  In  doing  so  we  would  ignore  the  contribution  (zero  or  not)  to  the 
vertical  deflection  (or  displacement)  due  to  other  parts  of  the  body  that 
we  neglect  by  integrating  between  specified  limits. 

The  correct  v.'ay  consists  in  writing  (TIMOSHENKO,  1970)  : 

w  »  I  Cz  .  dz  +  / (a) 

where  / (4)  is  a  function  of  4  only,  and  thus  a  constant  regarding  z, 
so  that  by  differentiating  we  obtain  again 


i  k 


By  choosing  an  appropriate  expression  for  /(a),  we  obtain  the  bottom 
fixed  at  the  desired  depth  :  by  doing  this,  we  introduce  in  fact  a 
geometrical  condition  fixing  the  reference  level  for  the  vertical  deflec¬ 
tions  at  the  chosen  depth. 

But  with  the  general  solutions  of  the  compatibility  equations  in  multi¬ 
layer  theory,  v/e  can  do  it  also  in  another  way  by  expressing  that  w  =  o 
at  the  desired  depth  and  determinating  the  corresponding  values  of  the 
parameters  A^,  C^, 

Since  there  are  thus  several  possibilities  to  express  a  same  boundary 
condition,  it  is  necessary  to  compare  the  results  obtained  and  retain  the 
one  that  seems  the  most  appropriate. 

To  do  this,  we  shall  consider  the  most  simple  case,  that  of  tlie  semi- 
infinite  body  ;  the  one-layer  case. 

Let  us  consider  a  semi -infinite  anisotropic  body  submitted  to  a  uniform 
distributed  vertical  pressure  at  its  surface. 

The  stressfunction  is 

^  pa  I  (A/""  -  +  re"""  -  dm 

'  O 

The  surface  boundary  conditions  (o2=p,  z=e)  are  induced  from  the 

relations  given  in  §  1.1.  for  the  stresses 

n  (1+n)  (Am^  H  Bm^]  +  u  (mn)  (Tim^  +  P6m^)  =  1 

n  (1+ii)  (Am-  -  Rm^ )  +  u-i  (ninl  (fim^  -  Pim^)  =  o 

Solving  this  system  for  R  and  P  ,  one  obtains  : 

n  (I-4)  (1+p)  Rm^  =  -  i  +  m(1+4)  (l+n)  Am^  +  ?ni  (»Hu)  Cim^ 

n  (I-4)  (kHii)  Pim^  =  1  -  ?n  (l  +  )i)  Am^  -  n  (1  +  j)  (n+u)  Tim^ 

The  next  step  depends  on  the  boundary  condition  that  fixes  the  deflection 
at  a  depth  fl. 


Referring  to  the  stress  ans  displacements  equations  given  in  paragraph  1.1 
the  condition  h>  =  o  at  a  depth  H  is  written 

n  (1+ji)  (A»rp  -  Bm  e  )  +  (CABre  -  Pim^e  )  =  0 

Replacing  B  and  P  by  their  values 

An,^  -  nlliw)  (Uy)  ^  2n4|H+ul,-^m//  j 

(1-4) 

.  p..„2  ,  U4(u+Li)^„4m//  2«4|«+p)  „-mtl  .  tt4(l+4)  |«+ii) , 

+  Cim  1-71-^,.,  ?  -  -fT^  ^  +  (^p)  (ra*-  ^  ’ 


4  (w-Hi)  1  -inill 
(1+p)  1-4  ^ 


4 

1-4 


-niH 

e. 


During  the  integration  proces,  the  value  of  m  tends  to  infinity,  so  that 
the  limit  expression  of  the  equation  becomes 

Iw  f  AmMl+y)^  e"'^  +  C4m^4  (m+p)  ^  1  =  0 

This  relation  can  only  be  satisfied  by  setting 
f  =  0  when  4  >  1 

A  =  0  when  4  <  1 


It  is  easily  shown  that,  in  the  case  of  an  isotropic  body,  one  must  always 
have  C  =  0,  because  of  the  factor  z  multiplying  f  in  the  stressfunction. 
Taking  4  <  1,  r  becomes 


n4(ii+p)  Cim^ 


-4m(/ 


l+p)2  -  (n+p)e 


0(1.  I  "mH  -4mf/ 

2|l+p)c  e 


(h+p)  II-4)  -  11+4)  l»7+p)e 


-  iimfl 


The  expression  of  the  deflection  at  the  surface  and  in  the  axis  of  the 
load  (Jo(mA)  =  1  for  a  =  0  1  is 

w  =  .  pn  [  —  ^  (  1  -  2h4(>hp)  Csm^  ]  dm 


• .  •  w  •  •  •  V*  4 
- J 


^2i$?^!_!29^59fD_b^_§D_§DGCCD!rl§^£_f*JD£^l9D 

The  vertical  deflection  at  a  depth  z  is  piven  by 

.  1+P  f  Jo (»w) . Ji (ma)  ,  •  2  wz  „  2  -mZi 

IV  =  .  pa  — — nt  —  I  n(l+p)  (^‘nrp  -  Rm^e  ) 

^  O 

.  Miln+iij^.o.  2  .Amz  n,  2  -imz,  ,  , 

+  ‘  ■  Him  P  )  1  dm 

We  now  choose  an  appropriate  function  /(^i)  so  that  iv  becomes  zero  for  z  =  If 

ti,\  1+h  ^  f  Jo  (om ) .  Ji  (mrt)  ,  ,1.  1  11  2  ml/  p  2  -nilf, 

/(''<)  =  -  .  pa  — ! — ' —  |m(1+ij)  (Am'^p  -  l,m^e  } 

o 

.  nilmii)^  /p,  2  imll  „  2  -im/l,  ,  . 

+  (C4m  e  -  V6in‘‘e  )]  am 

The  final  expression  for  the  deflection  is  then 

IV  =  lOj  +  /U) 

One  verifies  that  /(^)  is  indeed  only  a  function  of  k  and  that 
(V  =  0  for  z  =  II. 

For  the  reasons  developed  before,  one  of  the  parameters  A  or  C  must  be  zero. 
The  other  parameter  is  obtained  by  a  supplementary  boundary  condition 
(a  mechanical  condition)  ; 

-  =  0  at  the  depth  II 

-  a  =  0  at  the  depth  II 

If  we  still  suppose  4  <  1  and  thus  A  =  0,  one  obtains  in  the  case 
that  T.  =0 

KZ 

-llni  -6llm  -?6llm 
e.  .  c  -  c 

)'I4(m+ij)  C6m^  =  - rr: - -77- - 


The  deflection  at  the  surface  is  then 


=  pa 


Ml-n)  fji 

f(l-a)  J 


lamj 


m 


dm 


+  pa 


f  J|  (am)  f  a  («+li|  .-al/m  a  .-mil 

J  m  ^  (1-a)  (1+p)  ^  ■  1-a  ^ 


+  na  (nn.)  Cam  (  (i.],  ^ 


*'>+11  .  ^^dm  H+y  ,  1+s  ^-Anill  ,  •,  , 

-  ui7“  ■  ThT"  -FJ  '  >  > 


In  the  case  n  =  o 
niln+y)  Ca»i^  = 


c.  I  (M+ylae  -  (l+y)g  I 


2a(H+y)e''^'"  +  U+y)  (1-a)  -  (1+y)  (l+a)?'^^"'" 


The  deflection  at  the  surface  is  given  by  the  same  relation  as  above 
with  the  appropriate  value  for  C. 

One  sees  that  in  the  3  cases,  the  deflection  is  composed  of  a  first  term 


-  M  l-nl 

'  ril-a) 


am 


m 


dm 


=  pa 


a  1 1-w) 


This  term  is  the  deflection  on  top  of  a  semi-infinite  body. 
In  the  case  of  an  isotropic  body,  one  has  («  =  1): 


10- 


pa 


F 


The  second  term  lo^  depends  on  the  choosen  boundary  condition  ;  but  in  this 
3  cases  it  reduces  the  value  of  w-  because  of  the  fixed  bottom. 

We  have  computed  the  values  of  w-  and  for  different  values  of  ll/a. 

The  results  are  given  below,  at  a  factor  (l+y)/E  ,  for  a  =  0.5  and  for 
the  isotropic  case. 
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u 

o 

w-  =  1.025 

HIc. 

10^  (u  •  ol 
n 

w  (w  =  0 
n 

1 

0.339 

0.729 

n.i53 

2 

0.120 

0.295 

0.043 

3 

0.058 

0.147 

0.019 

4 

0.034 

0.086 

0.011 

5 

0.022 

0.057 

0.007 

6 

0.015 

0.040 

0.005 

7 

0.011 

0.029 

0.004 

8 

0.009 

0.023 

0.003 

9 

0.007 

0.018 

0.002 

10 

O.OOfi 

0.015 

0.002 

n  =  1 

w-  =  1.000 

1 

0.402 

0.934 

0.276 

2 

0.143 

0.442 

0.086 

3 

0.069 

0.231 

0.040 

4 

0.040 

0.138 

0.023 

5 

0.026 

0.091 

0.015 

6 

0.018 

0.064 

0.010 

7 

0.013 

0.048 

0.008 

8 

0.010 

0.037 

0.006 

9 

0.008 

0.029 

0.005 

10 

0.007 

0.024 

0.004 

We  conclude  that  from  a  depth  of  about  U/a  »  5  •  the  absolute  influence 
of  the  fixed  boundary  is  neglijjible.  This  influence  will  still  be  much 
lesser  in  the  case  of  a  roadstructure  where  the  E-moduli  of  the  layers  are 
sensitively  higher  than  the  modulus  of  the  siibground. 

We  conclude  also  that  the  relative  influence  on  the  deflection  is  much  more 
important  when  we  fix  the  horizontal  displacements  («  =  o) .  This  should  be 
the  case  in  a  laboratory  testpit  with  lateral  walls,  but  less  in  the  case 
of  a  real  road  where  lateral  movements  are  not  restricted. 

The  relative  influence  of  the  condtition  t..  *  o  is  also  more  important  than 
that  of  the  condition  w  -  o,  altough  less  important  than  the  condition  u  ^  o. 
It  seems  nevertheless  very  unlikely  that  there  would  be  no  friction  between 
the  subground  and  the  last  layer. 

The  easiest  way  to  fix  the  bottom  from  a  mathematical  point  of  view  is  on 
the  other  hand  the  condition  lo  »  o. 

Taking  then  into  account  the  little  influence  of  the  choosen  condition  on  the 
deflection  at  the  surface  and  the  fact  that  conditions  u  =  o  and  »  o 
have  less  physical  sense,  we  retain  the  condition  wi  «  o  as  the  most  indicated 
fixed  bottom  condition. 


/ 


1.3.  SOLUTION  OF  P/>RTICULAR  NUMFRICAL  PROBLEMS 


The  value  of  any  stress  or  displacement  is  obtained  from  one  of  the  above 
mentionned  relations. 

Let  us  consider,  for  example,  the  vertical  stress  in  the  i-th  layer  of  an 
isotropic  layer  : 


0^  »  pa  f  Jg  (nvi) .  Ji  (ma|  f -  0' m(  l-?|i( 
+  £!tm(  1  dm 


The  integration  can  only  be  performed  numerically. 

One  must  thus  calculate  the  value  of  the  stress  for  a  set  of  values  of  m 
growing  from  o  to  a  value  high  enough  to  ensure  convergency. 

For  each  value  of  m,  those  of  the  parameters  A/,  Ui,  Ci  and  Vi  must  be 
determined  out  of  the  set  of  boundary  conditions,  a  system  of  {4n  -  1) 
equations  with  (4n  -  1)  unknowns  in  the  case  of  a  fixed  bottom  and  n  layers 
above  it. 

The  first  programs  solved  this  problem  by  inverting  the  matrix  of  the 
{4n  -  1)  unknowns.  Nevertheless  the  inversion  procedure  leads  in  some 
cases  to  unsoluble  difficulties  because  of  the  presence  of  the  negative 
exponents  tending  to  zero  in  the  determinant  of  the  denominator. 

Other  programs  have  tried  to  avoid  the  inversion  procedure  as  follows  : 
one  chooses  appropriate  values  for  £?«  and  On  ,  goes  through  the  whole  set  of 
equations  and  verifies  in  how  far  the  surface  conditions  are  met. 

One  then  chooses  another  pair  of  values  for  and  and  follows  the  same 
procedure.  Since  the  whole  process  is  linear,  the  correct  values  for 
and  Vn  can  finally  be  obtained  by  linear  interpolation  after  two  runs. 

The  difficulty  lies  in  the  appropriate  choice  of  the  values  of  B^  and 
Vfi  to  ensure  a  numerically  correct  interpolation. 

However,  even  those  programs  are  not  entirely  appropriate  for  the  cases 
with  frictionless  conditions  at  some  interfaces. 

We  shall  show  this  with  the  most  simple  case,  that  of  a  two  layer. 

Writing  Ai,  Bi,  Ci,  Vi  in  stead  of  B^m*,  Qm,  Vim  ,  the  boundary 

conditions  are,  in  the  case  of  two  isotropic  layers,  with  the  origin  (z  •  o) 
at  the  interface,  the  thickness  of  the  first  layer  being  H  and  the  second 
layer  semi-infinite  : 


At  the  surface  (z  »  -H)  : 


+  Bie'^  -  Ci  (l-2ui+mH)  +  Vi  ll-2ui-niHU"’^  =  1 
Aie’"^  -  +  Cl  (2iJi-mHlc''^  +  Pi  (2ui+niH)c'^  -  o 

At  the  frictionless  interface  (z  -  o)  : 

Ai  +  8i  -  Ci(l-2uil  +  Pi(l-2uj)  *  82  +  P2(l-2u2) 

Ai  -  Bi  +  2uiCi  +  2yiPi  »  0 

-  82  +  2u2p2  ®  0 

[A,  -  81  -  2Ci(l-2yi)  -  2Pi(l-2yi)l  -  (  -  82  -  2p2  ( 1-2^2)  1 


Solving  the  system  for  Ci  and  Pj  ,  one  obtains  (BURMISTER,  1943) 

Cl  «  [  (1- F-^mh')  -  (1  -  Flc'”'^]  .  ^ 

Vi  *  (  Fe'^  -  (F- 


where 

V  •  Fe^'"^^  +  (2F-l).2mH  -  |l+2m2H2|  + 


p  ,  (l-y2)  +  n(l-yi) 
2(l-y2) 


n 


.li.Il+HiJ 

EJ  (I+U2) 


Ai  -  CilF-2yi)  -  Pi(l-F) 


8,  «  C,F  -  Pi(l-2yi-F) 


The  vertical  deflection  at  the  surface  is 


w  ■  pa 


m 

iljf  Jo  (w^)  -Ji 

•1  J  m 


('”al  r  t  0  2'^' 

-  (A,ni^c  -  Zitrre 


-  l2-4wi+mH)  -  (2-4uj-mHl  dm 


Replacing  Ai,  Bj,  Ci  and  P,  by  their  values,  the  deflection  becoires 


E,  J,  »' 


-  {2F-l-2mf/) 


:i-F)e 


Fp^  +  |2F-l)2mM  -  (i+2m"H^}  +  (l-F)e  " 


At  the  origin  of  the  integration  (m  »  o),  the  tenr  in  between  brackets 
becomes  indefinite  :  C/0. 

This  has  no  influence  when  computina  stresses,  because  the  Pessel  functions 
products  occuring  here  are  also  zero  at  the  origin  :  Jo  (fta)  . Jj  |ma)  ^  o 
for  m  »  0  . 

But  in  the  case  of  the  deflection  Jo .  Ji  (ma)  ^  a 

irr^o  m  Z 

It  is  therefore  absolutely  necessary  to  have  the  term  in  between  brackets 
in  a  sufficient  closed  form  to  be  able  to  determine  its  value  for  m  «  o. 

The  importance  of  the  first  term  of  the  series  is  not  neglioible  : 

for  m  »  0  the  term  in  betv/een  brackets  is  equal  to  Fill-ull/IEzd-u,*)  1  . 

If  li  is  the  interval  chosen  for  the  numerical  integration,  one  can  then 
wri  te 


i'.'  =  -  pa 


2(l-uid  f  1  F,(1-1J2)  h  .  Jo(rrvr).J,(ma)  .  .  .  1 

i  ^  rTFw  I  ^  *  1  - ' 


and,  if  v/e  make  a  semi-infinite  body  from  the  two  layer  (Fi  ^  F2,  ui  »  uz)  : 


2(l-u')  f  h 


Jo  (m>t)  .J,(ma)  ,  , 


Comparing  this  expression  with  that  for  the  deflection  at  the  surface  of 
a  semi -infinite  body 

tOm  ^  -  pa.2ll-ij^}  /  E 

One  concludes  that  the  contribution  of  the  first  term  h/6  is  indeed  not 
negligible,  especially  when  we  have  in  mind  that  the  only  practical  measu 
rement  that  can  be  performed  on  a  real  roadstructure  is  the  vertical 
deflection  at  the  surface. 


1.3.2.  Over  and  underflow  problems 


During  the  integration  procedure  m  varies  from  o  to  a  value  high  enough 
to  ensure  convergency.  We  mean  by  this  that  the  integration  procedure 
can  be  stopped  from  the  moment  on  that  the  terms  of  the  series  become 
so  small  that  they  have  no  more  influence  on  the  final  result  and  can 
thus  be  neglected. 

Practically,  however  this  means  that  m  can  reach  quite  hioh  values  such  as 
20  or  30  for  example. 

To  illustrate  the  influence  of  this,  let  us  go  back  to  the  two-layer 
developped  in  the  preceeding  paragraph. 

The  values  of  Ci  and  Pi,  from  which  the  values  of  all  the  other  parameters 
can  be  deduced,  are 


f  (l-F-nnH)e'^  -  d-Fle'"'"  1 _ 

+  l2F-l).2mH  -  (l+2m*H*)  + 


I  Fc"”  -  1 

Fc^*^  +  (2F-l).2mH  -  lU2m^H^)  +  (1-F)e‘^"'” 


The  geometrical  units  are  generally  expressed  in  function  of  a, 
the  radius  of  the  load. 

Let  us  consider  H/a  »  5. 

One  immediately  sees  that  no  computer  can  handle  exponents  as  and 

e  '  without  overflow  occuring  for  values  of  m  above  10. 

However  this  problem  can  easily  be  solved  by  dividing  both  numerator  and 

2iitW 

denominator  by  e  ; 

(  (l-F+mH)e''"^  -  (l-F)e"^"'”) 

‘  F  +  [  (2F-l).2mH  -  ( U2m^H^)1  +  (l-F)e''’’"^ 


_ [Fg-'^  -  1 

F  +  [  (2F-l).2mH  -  ( l+2m^H*)l  +  (1-F|e‘ 


Underflow  will  obviously  occur  now,  but  most  of  the  computers  have  a 
routine  that  sets  variables  subjected  to  underflow  equal  to  zero. 

If  such  a  routine  does  not  exist,  it  is  very  easy  to  build  it  into  the 
program. 

But  more  interesting  is  the  fact  that,  having  transformed  the  relations 
for  Cl  and  Pi,  convergency  will  occur  quite  quickly  and  in  a  complete 
safe  way  :  the  numerators  tend  both  to  zero,  while  the  denominator  tends 
to  a  constant  F, 


This  can  be  obtained  automatically  in  writina  the  boundary  conditions  at 
the  surface  (z  =  -H)  as  follows  : 

-  Cl  +  Pi  ( l-2vi,-mH) 

-  Bie’'"^  +  Cl  (2yi-mH)c'^"'^  +  Pi  (2ui+mH)  =  0 

However  this  is  only  true  in  the  case  of  a  two-layer  system. 

In  a  three-layer  with  Hi,  the  thickness  of  the  first  layer,  and  Hj.  the 
thickness  of  the  second  layer,  occur  exponents  such  as 


.ZmHi.ZwHz 


-2mHi.2n-.H2 


but  they  eliminate  when  writinq  the  denominator  in  close  forme  so  that 

P.tnH  2rnH 

dividing  the  expressions  by  the  largest  out  of  e'  'and  e  ^is  enough. 
If  one  should  divide  by  ^  the  denominator  v;ould  also  tend  to 

zero,  which  should  stop  the  program  because  of  dividing  by  zero. 

Here  is  another  reason  for  writina  the  equations  in  closeform  for  three- 
and  more-layers. 


'  •.**  •s."  ' 


Kjr.fi 


•  B  •  f  *  #  -  I 

^ A  , 


.•-V- 


‘ 

*  .  '  -  ^  • 

•  •  •  •  K  •  ‘ 

u  v_v. 


■  m 


r. 
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1.3.3.  The  vertical  deflection  at  the  surface 

When  one  computes  the  deflections  at  the  surface,  converqency  is  obtained 
only  very  slowly. 

To  illustrate  this,  let  us  look  at  the  expression  of  the  deflection  at  the 
surface  developped  in  §  1.3.1. 

w  •  -  pa  f  Jo  (iwaI  .  Ji  (mg) 

E,  J,  m 

(  -  (2F-l-2mH|  -  d-Fle'^"^  .  ^ 

Fe^*^  +  (2F-l)?mH  -  (l+2m^H^}  +  (1-F) 


To  avoid  overflow  problems  we  divide  numerator  and  denominator  by  ^ 

ui  •  -  pa  I  [  Jo  •  Ji  (ma| 

El  J,  m 

,  F  -  i2F-l-2mH}(i'^"^  - 
^  F  +  t  (2F-l)2mH  -  (l+2m'M2)le'^"’^  +  (l-F)e''’'^ 


For  large  values  of  m,  numerator  and  denominator  tend  both  fo  F, 

so  that  for,  let  us  say  m  »  ,  the  expression  above  could  be  written 

as  follows  : 

"'l 

M  •  -  pa  f  Jo  (mg) .  Ji  (ma) 

f  1  •'o  m 


F  -  (2F-l-2mW)e'^'"^  - 


F  +  f  (2F-l)2mH  -  (U2m^H^}]  ^  ^ 

-  pa  [  Jo  •  Ji  lmfl) 

Fi  J  m 


dm 


(1) 


The  first  integral  converges  fast,  the  second  converges  proportionnaly  to 
1/m.  This  means  that  if  one  should  want  a  result  correct  at  10'®,  one  has 
to  perform  the  numerical  integration  of  the  second  integral  until  values 
of  m  above  100,000  ! 


22. 


This  is  practically  impossible. 

One  could  of  course  tempted  to  interrupt  the  integration  procedure  when  the 
first  inteoral  has  converged  toa  satisfactory  level,  "hoping"  that  the 
second  integral  can  be  neglected  at  that  moment. 

To  illustrate  the  danger  of  such  an  approach,  we  return  to  the  semi¬ 
infinite  body. 

In  the  case  of  an  isotropic  body,  the  deflection  at  the  surface  and  the 
vertical  stress  at  a  depth  z  are  given  in  the  axis  of  the  load  by 


w  *  -  pa 


2(l-u*) 


Ji (ma) 


dm 


(2) 


1 


*  pa  I  Ji(ma]  (1-mz)  e  dm 


(3) 


Those  integrals  can  of  course  be  solved  analytically 

-  p.  2(1V > 


(4) 

(5) 


We  can  compare  (2)  with  the  second  integral  of  (1)  and  (3)  with  the  first 
integral  of  (1). 

We  perform  then  a  numerical  integration  of  (2)  and  (3)  and  stop  the  pro¬ 
cedure  when  (3)  has  converged  *0  a  satisfactory  level,  which  is  easily 
checked  by  comparing  the  obtained  result  with  the  correct  one  given  by  (5). 
The  difference  between  the  numerical  result  for  w  obtained  at  that  moment 
by  integration  of  (2)  v;ith  the  analytical  result  given  by  (4)  will  give  us 
an  illustration  of  the  possible  error  v/hen  integra'^ing  (1)  and  stopping  the 
process  when  its  first  integral  has  converged. 

This  difference  is  illustrated  on  figure  1. 

In  abciss,  v/e  have  the  convergency  level  adopted  for  the  vertical  stress  and 
in  ordinate  the  error,  expressed  in  on  the  values  of  the  vertical  stress 
and  the  vertical  deflection. 

One  sees  that,  for  even  such  lov/  levels  as  10’,  the  value  of  the  vertical 
stress  is  absolutely  correct,  while  the  erroron  the  deflection  varies  between 
-  5%  and  +  8%,  depending  on  the  chosen convergencylevel  and  the  relative 
depth  at  which  the  vertical  stress  is  computed;  the  worst  of  all  is  that 
we  have  no  means  to  predict  either  the  direction  either  the  amplitude  of  the 


error  on  w. 


The  only  way  to  solve  the  problem  satisfactory  is  to  split  the  expression 
of  the  deflection  in  another  way  than  the  ore  we  had  done. 

We  write  first  the  expression  of  the  deflection  with  neoative  exponents 
only  : 

m;  =  -  pa  r 


F  -  -  ( 1-F)  _ 

+  (  (2F-  1  )2mH  -  (l+2(r^H"ll  +  (l-F)c’^"''^ 


]  dir 


We  divide  then  the  numerator  of  the  term  in  between  brackets  by  the 
denominator  : 


VO  =  -  pa 


2(l-ui^l  f  Jo  (nvi)  .Ji  lira) 


1  +  2 


(  (1-F)  (1+mH)  +  1  e 


2i,2i  .-2mH  ,, 


:i-F)e 


-2n’H 


F  +  (  (2F-l)2mH  -  (l+2m2H^l]  e.  +  !1-F)e 


-4mH 


dm 


and  split  the  integral  into  two  parts  from  v/hich  the  first  is  inteqrable 
analytically  and  the  second  converges  in  the  usual  way. 


For  /I  =  0,  one  has 
For  a  =  a,  one  has 
For  K  <  a,  one  has 


Ji (mg) 

III 


dm  =  1 


=  ?/n 
m 

Jo  (m/tl  ,  f  (j/2,  .j/2;  1;  Jll 


where  F  is  the  hypergeometric  function  of  G4USS  : 
F  (a,  b;  e;  z)  =  I 


(a)..  =  a  (g+l)  (a+2  )  ...  (a+n-l| 


For  K  >  a,  one  has 


Jo  (m.l)  .Ji  (r^l  ^  «  F  (//2^  112;  2;  —  ] 

m  ^2 

The  obtained  result  will  now  be  correct,  while  converoency  is  reached  as 
fast  as  for  the  other  equations  for  stresses. 

But  again,  if  we  will  be  able  to  cotrpute  as  indicated,  we  must  dispose  over 
the  equations  in  closeform,  altcuqh  in  such  a  form  that  the  integral  can 
be  split. 


1.3.4.  Accuracy  In  the  first  layer 


As  for  the  deflection  at  the  surface,  the  numerical  computation  of  the 
stresses  in  the  first  layer,  in  fact  nearby  the  surface,  converges  also 
very  slowly. 

To  illustrate  this  let  us  look  at  the  relation  for  the  vertical  stress  in 
the  first  layer  (the  equation  is  given  in  §  1.1)  of  a  two  layer  : 

=  pa  [  f  +  Bim^e  -  Cim(l-2ui-mz)c*”^ 

'o 

+  Pim(l-2jji+mz)e  )  dm 

We  replace  Ai,  8i,  Ci,  Px  by  their  values  obtained  in  §  1.3.1. 

«• 

Oz  ■  pa  I  Jo  (ma) .  Ji  (ma)  f  F  I  1  +  m{H+z)  )  t. 

...  f  ..'"2  .  f  i  r  --mz  ]  dm 

+  [  je  .e  +[  .e  +I  le  .e  T 

The  values  of  z  are  negative  in  the  first  layer  (z  «  o  at  the  interface). 
The  term  F  ( 1  +  m(K+zl  1  e*  .e”  .  converges  very  slowly  for  values  of 

-  z  nearly  equal  to  H,  the  other  terms  converge  normally. 

To  solve  the  problem  created  by  the  first  term  we  divide  again  numerator 
by  denominator  : 

F  (  1  +  m(H+z)  1  (  1  +  m|W+z)  J 

[1  *  m(H^■z)]  [  [  (2F-l)2mH  -  (U2m^H^)  )  (1-F)e‘^"^  ] 

F  +  (  (2F-l)2mH  -  +  d-Fle'^"^ 

The  second  term  of  the  second  member  converges  normally  so  that  we  can 
again  split  the  integral  in  several  parts  from  which  the  one  that  converges 
slowly  is 

pa  f  Jo  (mA.) .  Ji  (ma)  {  1  +  m{H+z)  )  dm 

•  A 


This  integral  is  known  as  a  LIPSCHITZ-HANKEL  integral,  but  only  some  parti¬ 
cular  cases  are  integrable  analytically.  To  solve  the  problem  for  all 
cases  we  have  to  make  a  detour  through  the  analysis  of  stresses  and  dis¬ 
placements  in  a  semi-infinite  body  submitted  to  an  isolated  local  force  P. 
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^  ^  m^J,(m/L)e-'"^  rfm  =  -^  -lii^ 


"  - 


'*  2(1-W^)  [ 
^  E  j 


?Tr  (^^+z^)y2 


Jo(m^)e  dm  -  iliHi  ni  Jo(m4)e”'"^  dm 

o  ^  e 


P  2(l-tJ^)  1  _  P  (Uy) _ ^ _ 

^  E  U'+ziM'/i  ^  f  Iz’m’IV! 


^.P.  I.UtO  jl-Eiij  f  *,  .  EJ  Qllil  f  „  J,|w)e-"'  rf, 

271*"  J.  2tt  E  4 

_  (Uy)  (l-2y)  ^  (^^+z^}'/2  z  _  jP_  (Uy)  zK 
277  E  n  (/l*  +  Z^)*/2  277  E  (z*+  4*1^2 


1.3.6.  Stresses  and  displaceiT’ents  under  a  uniform  distributed  load 


The  stresses  and  displacements  under  a  uniform  distributed  load  can  be 
obtained  by  integrating  the  relations  under  an  isolated  load  over  the 
concerned  area 


•277  rCi 
Op  »  Oj  p  d  p  d 


where  Oj  is  given  by  one  of  the  relations  of  §  1.3.5.  wherein  the  distance 
x  must  be  replaced  by  {x^+p^-  2xpcoi  6  )  */2 


■  n"\*  V 

%  V  ^ 

] 


*  r  -* 


t  / .  .-  -J 

MPW 


.■'■  ••■  .’■  -■• 


•  •  •  V.'  ■  - 

V.-.rV.V 
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*)  B®l55i2G§_f2r-5ti?.§^!r?§§.2^ 


The  relation  for  o^  under  a  distributed  load  is  given  by 

r  f 

o  ■  pa  1  JoifM.)  •Jiiiivi]e.  dm  -  pa  Jo  (m^l  .Ji  (ma)mz  e  dm 

""  J. 

«•  • 

-  pa  (l-2u)  Ji(ma)c‘"’"  dm  +  pa  f  J,(ma)2  e''"^  dm 

J,  fWl  J.  A. 


so  that  comparing  with  the  expression  for  the  stress  under  an  isolated 
load  we  can  conclude  that  : 


f  .  T  o 

11  :  pa  Jo  (mA.) .  Ji  (male  dm  =  —  - 

J,  2ii  J.  (22+^2+ 

12  :  PO-  Jo  (ma) .  Ji  (ma)mz  e*”'^dm 

O 

.  j-  2zp _  2f 

i  i  L  (z^+a^+o^-2ap  C04  e  )  V2  (r*- 


p*-2ap  coa  6  ) 


d  p  d  0 


3z(a^+p*-2ap  coa  0  )p 
(r*+A.*+p*-2Ap  coa  0  )  ®/z 


d  p  d  6 


Ji  (ma)  e  dm  * 


2ir  ,a 


p  r 

1. 


(z^+a.^+p^-2ad  coa  0  )  ‘/2  -  j 


.  K  (a*+p*-2ap  cca  0  )  (z*+a*+p*-2a.p  coa  0  r  * 


p  d  p  d  0 


I,  ;  pa  f  2l^.j,{^a)z  c-'”"dm  . T 

^  2  J.  J.  (z^+A*+ 


(z^+A*+p*-2Ap  coa  0  )  Vz 


B§l§^l22§-f2!r_t!]2_I^C§§!_2z 

These  relations  can  be  deduced  from  those  established  for  the  stress  o  . 

A 


c)  Bgl3^i9^!_f2C-^!3§-i^C£§§-I■tz 


m 

•  pa  I  Ji  (nw.)  .Ji  (malm  z  dm 

Is  :  pa  f  J.(^).J.(m)mze'’"^dm  ■  £-  f"'  f 

2  4  J,  (2*+>i*+P*-2^p  C04  6  1*^2 


‘^)  5Sl55i2Dl.f9C-5t!S_^5tli£9l-^i5Bl59®?®Ci-‘i’ 

«,  -  -  pg  glllHl’  f 

E  J,  m 

m 

-  pa  l~il  I  Jo(ma)  .Ji  (mtt)  z  e**"^  dm 


I,  :  pa 


f  Jo (mg) .Ji (mft)  . 

J.  m 


-mz,  ,  Pdpd 


p  dp  d  6 


(Z*+a*+p*-2^P  C04  0  )  */2 


I7  ;  pa  f  Jo(ma) .  Jj  |ma)  ze*"'*dm  • 


z*Prfp. 

i-  i.  tr2+7,*4.r,2.? 


de 


K  K  (z*+a*+p*-2/tp  coaO  I  Vi 


®)  Bsl35i222_f22-522_!}2!!;i?2259l-2i5Bl22S!!!?25-2 

U  .  pa  f  ^-mz 

E  J,  m 

-  pa  I  Ji  (nw) .  Ji  |ma)  z  e*"’^  dm 

m 

I.  =  pa  f  c-”'rf. 

■'•  m 

P  [*  (  (z^-»/t^-*-p^-2apcoAe  )  Vz  -  z]  pdpde 

^  (/i*+p^-2ap C04  0  ) '/2  (z*+7i^+p2-?7ip Poi  e  IV2 


P  z(^*+P*-2apcoa  0)‘ 

cun  ■  ■  ■  ■  '  '■" — 

?■''  j  |z*+/i*+p*-2>ip  COA  0) 


J,  ;  pa  Ji  lM).Ji  Im)  z  t 


p  ri  p  d  0 
IV* 


*7 

■  i-l 


1.3.7.  Resolution  of  the  double  integrals 

The  integrals  of  §  1.3.6  are  most  easily  solved  in  transforming  the  variable 
0  and  p  by  setting  :  x  «  p  coa  0  y  ■  p  A/n  0  pcfpcfe  •  dx  du 


One  obtains 


r 

2’f  J . 


2z  (a^-x^)V2 _ _ 

a  (z*+x®+a.*-2xa.)  |z*+A*+a*“2xA)  Va 


that  can  easily  computed  numerically. 

®  2z(^*+x*-2/ix)  |a*-x*)‘/*  r  1 _ _ 


-a  {z*+x*+/i*-2x/l1  (z*+>i*+a*>2JW.)  V«  l■(;l^+a*-2xA+z*) 


(z*+x*+^*-2«) 


[  *  2z(a^-x^) 


... 


a  {z*+A*+tt*-2xa) ’/a 


P_  f  “  _ L 

Z’’  J-a  (z^+x* 


2z  (a^-x^)V^ _ _ 


(z^+x*+a.*-2^x)  (z*+A.*+a*-24x|V2 


I3,  I*  and  I5  are  particular  LIPSCHITZ-HANKEL  integrals  that  we  shall  solve 
in  next  paragraph. 

I  .  JL.  f  ^  tn  t^^-2x>n-fl^+z^|  Va  ■»  (a^-x^|Va^j^ 

J-a  (/i*-2xa+a*+z*)  Vz  -  {tt*-x*)Va 


1 7  can  be  deduced  from  I| 

I*  and  Ij  are  particular  LIPSCHITZ-HANKEL  integrals. 


Resolution  of  the  LIPSCHITZ-HANKF.L  Integrals 


The  solution  of  the  Lipschitz-Hankel  Inteprals  Is  qlven  by  WATSON  {I960) 


1^  r  (u*2v)  r  f 

r  12v+l) 


(  }!*h,  v+1;  - 


where 


w*  »  6*  +  c*  -  2  be  coa 


jj  ,  E*  [  Ji (m^) . Ji (wfl) 
^  J  fft 


J_  f  F  (1.  3/2;  2;  -  Siif)  ain*  ^di 
z*  2tt  J,  z» 


with  (WAYLANO,  1970) 


F  II,  3/2;  2;  ^  Tl  -  (1  + 

jZ  yjZ  L  yZ  J 


and  w*  ■  a*  +  V  -  2  oA  coa 


z  f "  1  r 

The  resultinaintegral  Ij  'p  —  1  -  — - -  a-cn  ^d 

^  J.  w*  L  (z»w»)V*J 


can  easily  be  solved  numerically. 


m 

.  p  Ji!.-L  r  F  13/2,  2;  2;  )  a^n*  ^d< 

z»  TT  i 


F  (3/2,  2;  2;  -  (1  +  — ) 


w'^.-Vz 


*  *.'  %*  “,**"*  'V**-  '  .  -  -  •  .  •  V  *.*  *-• 


When  stresses  and  displacements  are  computed  in  the  axis  of  the  load, 
one  has  : 


Ii  :  pa  [  jj(ma)e  dm  * — — 
K  a 


a^+z^-z 


aHz^ 


I2  :  pa  f  Ji(ma)mze  dm  *  — - '  •  P<^ 

i.  (a^+z^'lVj 


I3  :  pa  I  pa  I  J.(-ra)c'"’"  dm  - 


Ii,  :  pa  f  .  Ji  (ma)  z  e"'"^  dm 

«  A 


1  [  -mz 

pa  Ji(ma)mzp.  dm  • 


Is  :  pa  Ji  (rriA.)  .Ji  (ma)m  z  dm  «  0 


u  :  pa.  f  iiJS)  .  pa 

I  in  u 

'o 


«• 

I7  :  pa  I  Ji  |ma)z  dm  »  z  .  Ti 


Is  :  pa 


J. 


(ma)  -mz  , 

^ — !  C  dm  *  0 


1 9  ;  pa  Ji  (ma)  .Ji  (ma)  z  e  "“dm  «■  a 


To  be  applicable,  all  the  developpments  of  paraqraph  1.3.4.,  again  requi 
all  the  equations  to  be  available  in  closeform. 


» -  r.  jr'.tnrir^  i-«  '.-v  ^ 'T’.-rt'r 


PART  2  :  NUMERICAL  RESOLUTION 


2.1.  THE  MATHEMATICAL  ANALYSIS 

We  give  here  only  the  main  principles  of  the  analysis.  The  detailed  analy- 
sisses  are  given  in  appendix. 


2.1.1.  Full  or  partial  friction  at  the  interfaces 
Next  system  of  boundavy  conditions  has  to  be  solved 
-  At  the  surface  -  P 


interface 

ZX. 

•  ’’zl  •! 

*  J 

«,w 

At  the  bottom 


w  =  0 

A^  or  Cu  *  0 


This  is  a  system  of  16  equations  with  16  unknowns  (Aj,  8i,  .  C„,  P^). 

It  is  very  difficult  to  solve  this  system  in  a  complete  analytical  way, 
as  BURMISTER  did,  for  example  for  the  two  and  the  three  layered  structures 
Indeed  it  is  nearly  impossible  in  doing  so  to  avoid  mathematical  errors  in 
the  analysis  because  of  the  extension  taken  by  the  different  expressions 
in  the  elimination  procedure  of  the  unknov/ns. 

Thus  it  is  necessary  to  reduce  the  analysis  to  a  more  comprehensive  model 
but  in  such  a  way  that  all  the  numerical  problems  detailed  in  part  1  can 
still  be  solved  completely  accurately. 


W: 


i 


i"  •'.•■%»  ■  '  J.'  ..TH*' 


•  .  !•'.  i' J  -■.  f. 


-'.  -Wjr 

36. 


The  main  objective  is  to  obtain  an  expression  for  each  unknown  parameter 
consisting  in  a  numerator  containing  negative  exponents  only  and  a  denomi¬ 
nator  containing  a  constant  term  and  neoative  exponents. 

During  the  integration  procedure,  when  the  variable  tends  to  infinity, 
the  numerator  will  then  tend  to  zero  and  the  denominator  to  a  constant 
value. 

It  is  thus  absolutely  necessary  that  all  the  exponents  appear  in  closeform 
in  the  analysis.  The  factors  multiplying  the  exponents  may  then  be  expressed 
in  a  more  comprehensive  form. 

The  mathematical  analysis  can  then  be  resumed  as  follows. 


First_steB 


Replace  in  the  boundary  equations  of  the  third  interface  the  parameters  A<, 
and  Cl.  by  their  values  obtained  from  the  fixed  bottom  condition. 

Secgnd^steg 

Write  all  the  interface  conditions  in  matrixform. 

-  At  the  surface  MI  (A,  8,  C,  (I  C)^ 

where  MI  is  a  2  x  4  matrix 

-  At  the  first  interface  Mi(A|  Bj  Cj  Pi)^»  M2(A2  Bj  Ci  ^2)^ 
where  Mi  and  M2  are  4x4  matrices 


k-’ 

!:• 

i,-- 

t- 

■ 

f-z 

L.-' 

L  • 

i 


i 


-  At  the  second  interface  M3(A2  8i  C2  *  M^|A3  83  C3 
where  M3  and  M.,  are  4x4  matrices 

-  At  the  third  interface  MsiAs  83  C3  Vj)^  =  Mj(?4 

where  Ms  is  a  4  x  4  matrix 
and  Me  is  a  4  x  2  matrix. 

Ihird^^steg 

Invert  the  matrices  Mi,  M3  and  Ms 
The  system  becomes 
MI  (Ai  81  Cl  Pi)^  »  (J 

(Ai  8,  Cl  Vi)^  .  Mr!M2(A2  82  C2  P2)^ 

(A2  82  C2  P2)^  «  Ms’IMmIAs  83  C3  P3)^ 


(As  Bj  Ci  V3]  -  M;}M6(84  Vu) 


and  finally 

MI.M;^M2.M3^M...Mi^Ms(B^  Pul^  «  (I  0)^ 

The  product  MI.Mi*.M2.M3‘ JU.Ms'.M*  is  a  2  x  2  matrix,  so  that  we  can  write 


1  ail 

2  1 

/  i 

1  ■ 

\  ^2  1 

<^22  i 

1  n.  j 

1  0  1 

and  solve  the  resulting  system  : 

0-22 

84  «  - 

^11.^22  "  fl21.^tl2 

-a2  1 

=  - 

0-\\>0.22  -  <l2  1  2 

I 

Utilizing  the  different  matrix  equations,  we  express  then  the  other  parame 
ters  in  function  of  84  and  r4. 

To  obtain  the  exponents  in  closeform,  all  the  matrices  are  split  in  such  a 
way  that  the  exponents  can  be  put  out  of  the  brackets. 

For  example,  in  the  isotropic  case  : 

Ml*  *  -  ^  j  ^  j  (  e  ’^.Mn  +  c’^.Mi2  ] 

M2  ■  [  e^.M2i  +  c  ^.M22  1 

Ms*  “”T(T-uTr  ^  ^  1  +  C'^M3  2  1 

M4  -  [  e^.M4i  +  e‘^.M42  1 

Ms*  ■  -  ^  J ;-y3 )  ^  ^  ^-Msi  +  C^.M52  1 

Me  »  (  e  ^.Mji  +  e.  ^,^(s2  1 
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where  x  *  mhi 

y  *  m[h\  +  ^2) 
z  *  m(fii  +  hz  +  hi) 
u  ■  m{hi  +  hz  +  hj  +  2  hi,} 
h^  being  the  thickness  of  layer  i. 

The  terms  of  all  the  matrices  are  of  course  known  in  closeform,  so  that 
they  can  be  introduced  as  input  values  for  the  numerical  procedure. 

The  terms  of  the  resulting  products  are  not  expressed  in  closeform  :  it 
is  sufficient  to  know  which  coloms  or  rov/s  of  each  matrix  contain 
only  zeros.  Indeed  in  the  final  product  intermediate  matrices  disappear 
because  they  are  identic  zero.  Practically,  those  matrices  vanish,  which 
are  preceeded  by  positive  exponents. 

The  values  of  all  the  unknov/n  parameters  can  then  be  expressed  in  the 
desired  v/ay,  a  numerator  vMth  only  negative  exponents  and  a  denominator 
with  a  constant  term  and  negative  exponents,  excepted  for  the  parameters 
81  and  Vi  at  the  surface. 

Indeed  at  the  surface  the  values  of  the  parameters  Bi  and  Pi  contain  in  both 
numerator  and  denominator  a  constant  term,  followed  by  negative  exponents, 
so  that  satisfactory  convergency  cannot  be  obtained. 

The  parameters  81  and  Vi  are  then  expressed,  utilizing  the  surface  conditions, 
in  function  of  the  parameters  Ai  and  Ci. 

The  resulting  expression  for  a  given  stress  or  displacement  becomes  then 

m 

a  *  pa  j  Jo  |m/i)  .Ji  (ma)  [  K  +  /i(Ai)  +  /2(ri)  1  dm 

'  O 

where  K  is  a  constant  term  and  filAi)  and  fzICj)  are  functions  of  the  para¬ 
meters  Ai  and  Ci  v/hich  converge  normally. 

The  expression  of  the  stress  can  then  he  split  into  two  parts 

m 

Oi  *  pa  Jo  Im'i)  .Ji  (ma) .  Kdm 

*  O 

which  can  be  solved  analytically 

m 

02  »  pa  Jo (mA) . Ji (ma)  [  /i(Ai)  +  /2(Ci)  ]  dm 

which  can  be  solved  numerically  in  an  accurate  way. 


The  same  difficulty  arises  in  the  first  layer  near  the  surface. 
Although  all  the  parameters  converge,  again  the  parameters  8i  and  Pi 
converge  very  slowly  so  that  complete  accuracy  is  also  here  hard  to 
insure. 

Therefore,  again  utilizing  the  surface  conditions,  we  split  the  expres 
sions  for  stresses  and  displacements  into  two  parts. 

The  first  part,  known  as  a  LIPSCHITZ-HANKEL  integral  (paragraph  1.3.4) 
is  solved  analytically;  the  second  part  is  solved  numerically. 


The  procedure  described  in  previous 
The  boundary  conditions  are  : 

paragraph  cannot  be  applied 

here 

-  At  the  surface 

0^  «  P 

-  At  the  first  two  interfaces  : 

'^Z4  *  ^Z-i  +  1 

(1) 

^>IZ4  • 

(2) 

T  ■,  *  0 

KZA.  +  1 

(3) 

A.  4+1 

(4) 

“  At  the  third  interface  : 


o  .  =  o  ... 

ZA.  Z.<.  +  1 

^KZi  '  ^^Z.C  +  1 

(•'.  «  W.  .  , 

A.  4+1 


•  U  .  « 

4+1 


At  the  bottom  : 


w  »  0 

or  Cl,  •  0 


The  systems  of  equations  at  the  first  two  interfaces  cannot  be  expressed 
in  matrixform,  because  two  of  the  four  equations  are  homogeneous. 

The  mathematical  analysis  is  here  : 

First_ste9 

Replace  in  the  boundary  equations  of  the  third  interface  the  parameters 
Ah  and  Ch  by  their  values  obtained  from  the  fixed  bottom  condition. 

Second_steg 

Write  the  boundary  equations  at  the  third  interface  (friction)  in  matrix 
IA3  B3  C3  P3)^  »  Ms '.Ms  (84  Ph)^ 


Ihird_steB 

Using  the  surface  conditions,  express  Aj  and  Pi  in  function  of  Ci  and  Pi 
(A,  8,)^  -  MolCi  P,)^ 

Replace  in  conditions  (1)  and  (2)  at  the  first  interface  Aj  and  8i  by  their 
values  and  solve  the  system  by  expressing  Ci  and  Pi  in  function  of  A2,  82, 
C2,  P2. 

(Cl  P,)^  .  M,(A2  82  C2  P2)^ 

Replace  in  condition  (4)  at  the  first  interface  Ai  Bi  Ci  and  Pi  by  their 
values  in  function  of  A2  82  C2  and.P2. 

Using  then  conditions  (3)  and  (4)  at  the  first  interface,  express  A2  and  82 
in  function  of  C2  and  P2 

IA2  82)^  «  M2(C2  P2)^ 

Replace  in  condition  (1)  and  (2)  at  the  second  interface  A2  end  82  by 
their  values  and  solve  the  system  by  expressing  C2  and  P2  in  function  of 
A3,  83,  C3,  P3. 

IC2  p2)^  •  M3(A3  83  C3  p3)^ 

Replace  in  condition  (4)  at  the  second  interface  A2  82  C2  P2  by  their 
values  in  function  of  A3  83  Cs  P3 

Equations  (3)  and  (4)  at  the  second  interface  reduce  then  to  following 
system 

M^(A3  83  C3  P3)^  =  (K  0)^ 

where  Mi,  is  a  2  x  4  matrix  and  K  a  function  of  the  integration  variable. 
One  obtains  finally 

M4. Mi'. We (84  Pu)^  *  (K  0)^ 

a  system  v/hich  can  be  solved  in  the  same  way  as  in  previous  paragraph. 

The  same  procedure  is  also  utilized  to  express  all  the  unknov/n  parameters 
in  function  of  8^  and  Pu. 
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Nevertheless  a  supplementary  difficulty  arises  in  the  full  slip  case, 
regarding  the  expression  of  the  vertical  deflection  : 


f' 


Jo  (iM) .  Ji  (ma) 


f.lA.  8.  f.  V.)  dm 

^  4,  4, 


When  integrating  numerically,  this  expression  becomes,  for  m  «  o 

tin  Jo  (ffVt) .  Ji  (ma)  ^  ^ 
m«o  m  2 


tun 

mso 


/.|A.  B.  C.  V.) 

X.  A.  A.  4.  A. 


0_ 

0 


SO  that  at  the  origin  (m=o),  the  expression  for  the  vertical  deflect! 
is  undefined. 

But,  thanks  of  the  fixed  bottom  condition,  we  know  that  at  the  bottom 
m;=o  for  all  values  of  the  integrating  parameter  thus  also  for  mao. 
Using  the  boundary  conditions  we  can  then  write  : 

At  the  bottom  /m(Am  8^  Cu  Pi.)  »  o 

At  the  third  interface  /3IA3  63  C3  Ps)  »  /i,(Ai.  Bm  C4  Pm) 

At  the  second  interface  fz[A2  82  C2  Pa)  =  /3IA3  B3  C3  P3) 
and,  at  the  first  interface  /i(Ai  81  Cj  PJ  =  /2(A2  82  C2  P2) 

Thus,  in  all  layers,  v/e  have,  for  m=o  :  /.(A.B.C.P.)  *0 

4,  ^ 
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2.1.3.  The  detailed  mathematical  analysis 
The  completely  detailed  analysis  is  piven 

in  appendix  1  for  the  isotropic  case  with  full  and  partial  friction 

in  appendix  2  for  the  isotropic  case  with  full  slip  at  the  first 

two  interfaces 

in  appendix  3  for  the  anisotropic  case  with  full  and  partial  friction 

in  appendix  4  for  the  anisotropic  case  with  full  slip  at  the  first 

two  interfaces. 


The  stresses  and  displacements  are  obtained  by  numerical  integration 

of  expressions  such  as 

•» 

CT  =  pa  [  Jo (m^) • Ji (ma)  /  (m  )  dm 


Therefore  we  use  Simpson's  method. 

I  la,b)  =  [/o  +  4/i  +  2/2 +. . .  +4/2M-1  +  fin  ) 

where  (a, 6]  is  the  interval  of  integration  of  I  (a,t)  subdivised  in 
2n  equal  segments  of  length  h. 

Here  the  interval  of  integration  goes  until  infinity. 

The  numerical  computation  is  interrupted  when  the  value  of  the  function 
f  lm,z)  (in  fact  the  values  of  all  the  parameters  A.  8.  C.  V .)  becomes  smaller 
than  the  imponed  convergency  level  : 

/^(m,zl<  e 

The  function  /^(mz)  varies  quickly  for  small  values  of  the  integrating 
parameter  m  and  much  slower  for  hinh  va;lues  of  m. 

Thus  it  seems  appropriate  to  increase  the  value  of  h  for  higher  values 
of  m. 

The  numerical  integral  can  then  be  written  : 

hi 

I  [o,»l  =  —  (/o+  A/i  +  2/2  +  4/3  +  /^  ]  for  m  <  Li 

3 

hi 

+  —  [  /^  +  4/5  +  2/5  +  4/7  +  fs  1  for  Li  <  m  <  Li 

3 

b-3 

+  —  [/b  +  4/9  +  2/jo+  4/ji+  /ij]  for  L2  m 

3 

+  ... 


For  practical  reasons  we  take 


Unless  in  the  case  of  the  vertical  deflection,  for  which  the  function  is 
computed  in  a  separated  v/ay,  the  first  term  fo  =  c  ;  the  last  one  f2n 
(here  fiz)  can  be  neglected. 


The  integral  becomes  then  : 


I  -y  [4/1+2/2+4/3+2/41+  ■5/4 


2h 


?.h 


+  "Y"  [^■[5  +  2/5  +  4/7  +  2/8  1  +  f 8 


4fi 


+  [  4/9  +  2/10  +  4/11  +  2/12  1 


The  main  computation  routine  is 


^  , 


I  U^,  1  ■  Y  '  *  2/^.,  *  ...  1 


When  modifying  the  leng^^of  the  integration  segment,  a  secondary  routine 
increases  the  already  obtained  with  a  value 

1”  ‘ 

The  initial  length  h  of  the  integration  segment  is  chosen  by  the  user. 
It  is  multiplied  by  a  factor  2  when  the  values  of  flmzj  become  smaller 
than  10’^  10■^  lO""  and  10'=. 

The  final  convergency  level  is  10'®. 


2.3.  RESUMED  FLOW  SHFET 


The  main  steps  of  the  program  can  be  resumed  as  follows. 

2.3.1.  INPUT  procedures 

Input  of  the  data  (by  display  or  file) 
regarding  loads  and  structure 

Input  of  the  coordinates  (by  display  or  file) 

of  the  points  where  stresses  and  displacements  are  to  be  computed. 
-  Choice  of  the  length  of  the  integration  segment. 


2.3.2.  Computation 

Computation  of  the  value  of  the  vertical  deflection  at  the  surface 
and  in  the  first  layer,  for  m  -  o. 

Computation  of  the  parameters  for  each  value  of  m. 

Computation  of  the  values  of  the  Pessel  functions  for  each 
value  of  m. 

Computation  of  the  stresses  and  displacements  in  cylindrical 
coordinates  for  each  value  of  m. 

Computation  of  the  stresses  and  displacements  in  cartesian  coordi¬ 
nates  for  each  value  of  m. 

Upbuilding  of  the  vectors  containino  the  results  using 
Simpson's  rule. 

Convergency  test  for  each  function. 

When  the  final  convergency  is  reached,  computation  of  the  additional 
values  at  the  surface  and  in  the  first  layer  (LIPSCHITZ-HANKEL 
Integrals) . 


OUTPUT  Procedures 


Creating  of  a  file  containing  the  results 
stresses  ox  Oy  Tx^  txz  "^yz 

displacements  ux  102 

Computation  of 

principal  stresses  oi  02  03 
principal  strains  Cj  Cj  ^3 
linear  strains  Cj,  Cy 

Printing  (if  desired)  of  the  results. 


Capabilities  of  the  program 


The  number  of  circular  loads,  with  different  radii  and  contact  pressures, 
is  limited  to  20. 

Stresses  and  displacements  can  be  computed  at  3C  places  in  the  horizontal 
plane  and  at  each  place  at  30  depths  (included  8  values  at  the  surface 
and  the  interfaces). 


2.3.5.  Used  language 


The  language  in  which  the  program  is  v/ritten  is  FORTRAN  77. 
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2.4.  PARTICULAR  FEATURES 


2.4.1.  Verification  of  the  accuracy  of  the  results 


Stresses  and  displacements  are  computed  at  each  Interface  v/ith  the 
appropriate  stress  function.  This  means  that,  at  each  Interface, 
the  computations  are  performed  at  the  bottom  of  layer  i,  usino  parameters 
8^  Ql  ,  and,  at  the  surface  of  layer  -c+1,  using  parameters 


1 »  1 »  J »  1 

Thus  stresses  and  displacements  are  computed  at  a  same  point  with  two 
different  functions.  This  procedure  allows  to  verify  the  results 
(mainly  to  insure  that  the  chosen  length  of  the  integration  segment  is 
small  enough)  takino  into  account  that  the  boundary  conditions  have  to  be 

satisfied  : 


-  In  the  friction  case  : 


02^  '  ^z-c*  1 

1 


-  In  the  full  slip  case  :  »  °zi*1 

This  is  of  particular  importance  at  the  first  interface  where  the  stresses 
at  the  bottom  of  the  first  layer  are  computed  in  a  complete  different 
way  (part  of  it  analytically). 

In  the  full  slip  case,  one  must  be  carefull  by  choosing  the  length  of 
the  integration  segment. 

Indeed  by  choosing  it  to  small  (0,02  for  example  when  a  good  average  value 
is  0,1)  there  is  a  risk  that  the  values  of  the  function  /  (mz) ,  undefined 
for  m  *  0,  result  in  complete  abnormal  figures  for  values  of  m  near  to  o  ; 
for  example  tensile  stresses  for  o^. 

Beside  the  proceeding  procedure,  the  best  way  of  verifying  the  results 
Is  to  compare  them  with  existing  tables  such  as  these  of . JONES  (1962) 
established  for  a  three  layered  system  with  full  friction  at  the  interfaces. 
It  Is  very  easy  to  reduce  the  four  layered  system  to  a  three  layered  one 
by  setting  the  Young's  moduli  and  Poisson's  ratio  of  two  successive  layers 
Identic.  That  creates  no  numerical  difficulties. 
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On  the  other  hand  it  is  necessary,  because  of  the  fixed  bottom,  to  oive 
an  important  thickness  to  the  last  layer  in  order  to  assimilate  it  to 
a  semi-infinite  body,  as  considered  by  Jones. 

Of  course  the  comparison  of  the  results  can  only  be  done  at  the  two 
interfaces,  but  there  also  the  values  of  Jones  are  completely  accurate 


2.4.2.  Converqency  of  the  Eessel  functions 


The  Bessel  functions  are  computed  by  series  expansion  such  as 
Jo(x) 


2  - 


,fe  (x/2)' 


fc :  h 


For  high  values  of  x  this  expansion  leads  to  overflow  problems,  because  of 
the  factorials  in  the  denominator. 

For  values  of  x  higher  than  16,  the  Bessel  functions  are  approximated  by 
their  asymptotic  values  given  by 


Jo(x) 

~  Va 

^  nx 

C06 

(X  - 

V4) 

Ji(x) 

=  yi 

7TX 

6-in 

(X  - 

^/4) 

“4 


The  introduced  error  is  less  than  10 
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2.4.3.  Underflov/  problens 

FORTRAN  77  does  not  have  a  routine  that  sets  variables  subjected  to  under¬ 
flow  equal  to  zero.  Thus,  this  routine  had  to  be  built  in  the  program. 

The  limit  value  for  underflow  on  the  IBf'  PC  is  exp  (-88) 

Maximum  there  occur  in  the  computation  4  products  of  exponentials. 

So  in  limiting  the  value  of  each  exponent  to  (-20)  normally  there  should 
occur  no  underflow. 

We  have  that  exp  (-20)=  lO'® 

Normally  one  computation  requires  about  60  loops,  let  us  thus  say  a  maximum 
of  100. 

The  maximal  error  resulting  in  the  limitation  of  the  values  of  negative 
exponents  is  then 

100  X  10“®  =  10"® 

which  is  the  value  of  the  convergency  level  of  the  whole  procedure. 

In  doing  so,  we  insure  thus  accuracy  of  the  results  at  a  level  of  10"®. 


'.■VU 
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2.4.4.  The  coinputation  of  the  principal  stresses 

The  principal  stresses  are  solution  of  followina  relations  (TIMOSHENKO 
and  GOODIER,  1961) 

+  -C2.5  -  -<’.3  =  0 

where  the  stress  invariants  are 
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^1  = 

a  +  a  +  a 

X  y  z 

X2  ' 

a  .a  +  0^.0, 

X  y  X  z 

^3  - 

o„.o  .a,  -  2t 

X  y  z 

0..a..  -  -  t2. 


'X"^y'^z  ^'xtj'  '^yz'  '^xz  ■  ‘  "^z'^^xy 


The  equation  is  transfomed  into 
+  of  +  6  =  0 

with  .t  =  S  +  .C1/3 

a  »  ‘^(3x2  -  xi) 
b  =  (-2'Ci^+9.ci.X2  ~  27x3) 

Knowing  the  trigonometric  identity 
4coa^e  -  3coa  0  -  coa(36)  e  0 
we  may  write,  v/ith 
t  *m.co6  6 

m^coa^0  +  amcoa0  +  fa  H4coa^0  -3coa0  -  coa(30) 
=  —  and  m  =  (-  4^)  *^2 


so  that  coa(30) 


The  solution  is  then  given  by 


e.  .i 

arc.  coal-^l 
a 

•  coa 

6i 

Hi  - 

-  4i/3 

e..i 

,3fa,  ^  2v 
arc .  coa  l~|  +  rr 

tt  3 

•  coa 

e* 

S2  ti 

-  -<1/3 

e..i 

i36,  .  4ir 
arc.  coal—)  +  — 

•  coa 

01 

53  *  ti 

-  ^1/3 

When  tt  -  0,  the  value  of  coalse)  seems  to  become  infinite. 
Nevertheless,  this  is  not  true. 

We  have  always 

K  -  V’”  ^ 

0 

l<’l  -  0,1  ’  >  0 


and  thus 


^  Sf  ^  V 


+  o. 


By  adding  o*+o*+a^  >  o  .a  +a  o  +  o  o 
X  y  z  X  y  y  z  z’  K 


Expanding  parameter  a,  v/e  have 


•  “’x  -  “‘u  -  -Zo„.o  -Zo,.o,| 


*  </  yz  z  X 


4  Io,.o„  +  o„.a,  +  a..o„  -  3t* 


-3t* 


-3t* 


Since 


we  have  that  alv/ays  a  <  o 


The  parameter  a  can  only  be  equal  to  zero  when 


*  V  * 


This  is  only  possible  when  ° 


T  «  T  =  T  *  0 
xw  t/2  ZX 


Thus  when  o^,  o^  and  are  the  principal  stresses  in  a  spherical  stress 
situation. 


The  value  of  b  is  then 


b  .  i  (  54  o*  -  27  o*  -  27  oM  =  o 
27 


and  the  value  of  m  is  undefined. 


Me  conclude  that,  when  a^,  the  principal  stresses  are  immediately 


obtained  by 


Oi  ■  02  ■  O3  *  '  o^ 


■  •.* 


r-.-.v'C-.-.-. 


^  •  In  •  • 


■ .  •  • '  -  ■  » 


??> 


a?- 


•  •  •  -.0  -L. 


►*•  ..■*• . 


fi .  - « • ./  . 


um. 


.'■  -■• 
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2.5.  USE  OF  THE  PROGRAMS 


2.5.1.  Generalities 

The  programs  are  available  in  two  versions  : 

-  EXECUTABLE  version.  In  which  the  executable  prooram  Is  made  up  of  only 
one  block  with  automatic  loading,  v/hen  DOS  Is  ready. 

-  SOURCE  version.  In  which  all  controls,  data  and  modules  are  in  separated 
files. 

Each  program  Is  v/ritten  on  a  separate  floppy  disk. 


The  names  of  the  programs  are  : 


FLIP  :  Four  layered  Isotropic  system  with  partial  friction 
at  the  interfaces. 

FLIS  :  Four  layered  isotropic  system  v/1th  full  slip 
at  the  first  tv/o  interfaces. 

FLAP  :  Four  layered  anisotropic  system  v/ith  partial  friction 
at  the  Interfaces. 

FLAS  :  Four  layered  anisotropic  system  with  full  slip 
at  the  first  two  interfaces. 


Each  diskette  contains  an  explanatory  notice  which  can  be  called  by 
next  instructions  : 


FLIP  NO. TXT  for  FLIP 
FLIS  NO. TXT  for  FLIS 
FLAP  NO. TXT  for  FLAP 
FLAS  NO. TXT  for  FLAS 


The  notice  FLIP  NO. TXT  is  given  In  appendix  5. 


V 
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2.5.2.  Source  version 


For  this  detailed  analysis  of  the  source  version,  we  refer  to  the  notice 
of  program  FLIP  in  appendix  5. 

Each  program  is  built  up  in  9  modules. 

Module  1  (FLIP  l.FOR) 

Main  module  controlling  the  whole  computation. 

Module  2  (FLIP  2. FOR) 

Data  input  module  with  following  subroutines  : 

DOMEC  :  Mechanical  data  of  the  structure 

DOTRA  :  Traffic  (load)  data 

POCAL  :  Stress  coordinates  in  the  xf/  plane 

POINT  :  Depth  coordinates 

PAS  :  Length  of  the  integration  segment 

CHECH  ;  Desired  geometric  scale 

ERROR  :  Error  routines  on  input  procedures. 

Module  3  (FLIP  3. FOR) 

Initialization  module 

ECHDE  :  Scaling  of  the  geometric  parameters 

VINIT  :  Initialization  of  the  integration  parameters 

ZERO  :  Initialization  of  the  results  vectors 

FINIT  :  Computation  of  the  values  of  the  vertical  deflections  at  the 
surface  and  in  the  first  layer  for  m  »  o 

Module  4  (FLIP  4. FOR) 

Computation  of  the  parameters  for  each  value  of  m 

P4442  :  Product  of  a  matrix  (4,4)  with  a  matrix  (4,2) 

PCT22  :  Product  of  a  constant  with  a  matrix  (?,2) 

PCT 42  :  Product  of  a  constant  with  a  matrix  (4,2) 

SOM 42  :  Sum  of  two  matrices  (4,2) 

SOM 22  :  Sum  of  two  matrices  (2,2) 

P4444  :  Product  of  two  matrices  (4,4) 

P2442  :  Product  of  a  matrix  (2,4)  with  a  matrix  (4,2) 

CONST  :  Computation  of  the  parameters  A-  C ■  B-  V- 


Module  5  (FLIP 5. FOR) 

Computation  of  the  stresses  for  each  value  of  m 

BESJl  :  Computation  of  Jj(tna) 

BJ0J2  :  Computation  of  Jo(m^)  and  Jiiom.) 

SURFA  :  Computation  of  the  stresses  at  the  surface 

COUC  1  :  Computation  of  the  stresses  in  the  first  layer 

C0UC2  :  Computation  of  the  stresses  in  the  second  layer 

COUC 3  :  Computation  of  the  stresses  in  the  third  layer 

COUC  4:  Computation  of  the  stresses  in  the  fourth  layer. 

Module  6  (FLIP  6. FOR) 

Integration  procedure 
TITRE  :  Display  control 

MODIF  :  Alteration  of  the  length  of  the  integration  segment 

SIMPS  :  Routine  for  Simpson's  rule  and  up  building  of  the  results  vectors. 

Module  7  (FLIP  7. FOR) 

Definite  values  of  stresses  and  displacements 

COMSU  :  Computation  of  the  additionnal  analytical  values  at  the  surface 

FONC  :  Computation  of  the  Gauss  function  for  the  vertical  deflection 

at  the  surface 

CCOU  1  :  Computation  of  the  additionnal  analytical  values  in  the  first  layer 
FONC  1  :  Computation  of  the  LIPSCHITZ-HANKEL  integral  Ij  (appendix  1) 

FONC 2:  Computation  of  the  LIPSCHITZ-HANKEL  integral  I2 

FONC 3  :  Computation  of  the  LIPSCHITZ-HANKEL  inteoral  T3 

FONC 4  :  Computation  of  the  LIPSCHITZ-HANKEL  integral 

FONC 5  :  Computation  of  the  LIPSCHITZ-HANKEL  integral  Is 

FONC 6:  Computation  of  the  LIPSCHITZ-HANKEL  integral  le 

ECHEF  :  Rescaling  of  the  displacements. 

Module  8  (FLIP 8. FOR) 

Printing  of  the  results 

IMOON  :  File  storing  and  printing  of  the  data 
IMRES  :  File  Storing  and  printing  of  the  results. 

Module  9  (FLIP  9. FOR) 

Display  instructions. 


^  nr 


2.5.3.  Running  of  the  prooram 


We  refer  again  to  the  notice  in  appendix  5,  on  paqe  5,  from  v/here  on 
the  complete  running  procedure  is  explained. 


2.5.4.  Miscellaneous 

It  is  important  to  notice  that  the  mechanical  input  data  (loadpressures  and 
Young's  moduli)  must  be  expressed  vvith  the  same  units  and  that  also  the 
geometrical  data  (radii  of  the  loads,  thickness  of  the  layers,  coordinates 
of  the  different  locations)  must  be  expressed  with  the  same  units. 
Nevertheless  it  is  not  necessary  that  mechanical  data  and  geometrical  data 
are  expressed  in  coherent  units  :  for  example,  mechanical  data  may  be 
expressed  in  psi  and  geometrical  data  in  meters. 

The  stresses  will  be  expressed  with  the  same  units  as  the  mechanical  data 
and  the  displacements  with  the  same  units  as  the  geometrical  data. 

-  Accuracy_gf_anisgtrogic_resylts 

It  is  not  possible  to  compare  the  results  of  the  anisotropic  programs 
with  those  of  the  isotropic  ones  in  setting  the  degree  n  of  anisotropy 
equal  to  one.  Indeed,  in  doing  so,  all  the  relations  become  undefinite. 
Approximate  checking  is  possible  by  setting  for  the  stresses  r.  =  1.05 
and  for  the  displacements  n  =  1.20  . 

-  Scale 

Routine  ECHDE  of  module  3  scales  the  Geometric  parameters  of  the  structure. 
This  routine  has  been  introduced  to  reduce  the  problems  regarding  accuracy  : 
in  using  the  same  scale  one  can  use  the  same  length  for  the  integration 
segement. 

Convergency  problems  arise  essentialy  at  the  surface  and  in  the  first  layer. 
The  best  scale  factor  is  therefore  the  thickness  of  the  first  layer. 

Another  appropriate  scale  factor  is  the  radius  of  the  loads,  obviously 
when  all  loads  have  the  same  radius.  The  computation  of  the  Bessel  function 
Ji(ma)  is  then  the  same  for  all  loads.  But  this  scale  factor  should  only 
be  utilized  when  its  order  of  magnitude  is  the  same  as  that  of  the  tickness 
of  the  first  layer. 

Of  course,  at  the  end  of  the  v/hole  procedure,  the  values  of  the  displacements 
have  to  be  rescaled  in  their  original  values  (routine  ECHEF  of  module  7). 
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Poisson's  ratio 


The  Young's  moduli  of  the  materials  are  knov/n.  This  is  not  always  the 
case  with  the  Poisson's  ratio.  Therefore  the  program  aks  if  Poisson's 
ratio  is  known.  If  not,  a  value  of  0.5  is  taken  in  the  isotropic  case. 
In  the  anisotropic  case  the  procedure  is  as  follows  : 

-  First  question  :  "Co  you  know  the  anisotropic  Poisson's  ratio  ?  " 

IF  YES  :  Input  the  value 
IF  NO  :  Second  question. 

-  Second  question  :  "Do  you  know  the  Poisson’s  ratio  when  considering 

the  material  as  isotropic  ?  " 

IF  NO  :  The  program  utilizes  the  maximum  tolerated 
value  :  0.5  for  an  isotropic  material 

for  an  anisotropic  material 

(LEKHNITSKII,  1963). 


IF  YES  :  Input  the  isotropic  value 


The  program  computes  then  an  adapted  anisotropic 


value 


,  l-5xn  ^ 
0.5  2  +  ff 
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The  program  computes  the  principal  stresses  in  routine  li’RES  of  module  8. 

It  also  computes,  in  the  isotropic  case,  the  principal  strains  by  applying 
HOOKE'S  law. 

It  does  not  compute  the  principal  strains  in  the  anisotropic  case. 

There  indeed,  the  values  of  the  Young  moduli  varie  v/ith  the  considered 
direction.  Thus  to  apply  HOOKE's  law  one  must  know  the  principal  directions 
to  be  able  to  compute  the  values  of  the  Young  moduli.  This  computation 
requires  the  resolution  of  a  system  of  three  trigonometric  equations  with 
three  unknowns.  Although  this  could  theoretically  be  done,  accuracy  of 
th  results  can  practically  not  be  insured. 


During  the  computation,  the  successive  values  of  the  integrating  parameter 
m  are  displayed.  This  is  only  done  to  show  that  the  program  is  running 
and  to  protect  the  user  from  despair  v/hen  utilizing  the  model  with  an  impor¬ 
tant  number  of  loads. 


"remain  patient"  on  the  display  means  that  the  propram  is  computing 
the  parameters  and  the  stresses  (modules  4,  5  and  6). 

"remain  more  patient"  means  that  the  program  is  computing  the 
LIPSCHITZ-HANKEL  integrals  (module  7),  which  can  need  a  certain 
amount  of  time  with  several  lo  ds  out  of  the  vertical  axis  where 
stresses  are  computed. 

2.5.5.  Values  for  n  and  1  . 

d 
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Granular  materials  and  soils  can  be  considered  as  anisotropic  due  to 
compaction  or  settlement.  In  most  of  the  cases  the  value  of  n,  the  degree 
of  anisotropy,  is  higher  than  1,  which  leads  to  the  phenomenon  of  vertical 
stress  concentration;  the  value  of  n  can  be  related  to  Frohlich's  stress 
concentration  factor  k,  approximatively  by 

2k  -  1  -\/4k  -  3 

n  a  - ^ - 

2 

The  value  of  k  varies  from  2  to  5  (3  =  isotropic),  thus  the  value  of  n 
from  0.4  to  2.4. 

For  superconsol i dated  clays,  n  is  lower  than  1. 

-  The  partial  friction  ratio  1.. 

Full  friction  is  introduced  by  setting  1^  =  1. 

When  1,  tends  to  infinity,  the  shearstresses  tend  to  zero. 

Thus  partial  partial  friction  is  obtained  for  values  of  1^  between  1  and 
infinity. 

The  practical  value  depends  on  the  value  of  the  maximum  shearstress 
(located  near  the  edge  of  the  load)  tolerated  at  the  interface. 
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APPENDIX  1 


Algebraical  Analysis  of  a  four-layer  isotropic  System  with  fixed 
bottom  and  partial  friction  interface  conditions. 


1.  Boundary  conditions. 


We  write 

Ai  =  Aiw’'  Bi  =  Ibino^ 

Ci 

Cw  -  X 

Ei(>t4^i) 

Fo.  - 

kw- 

Ww 

l_w  .  BltltL) 

Lw 

X  s.  Yvi  VI, 

y  «  IVj  W ,  4  ^  2-^ 

z  «  yv,  G  ^ 

t-  =  vv>  G»i-^^*‘-  ^ 

Hz) 

Iv  =  -  Z.V+  ^ 

=  C.rv,  X>|  - 
E-i  (-H  f*-*) 


£a(/l4')x3) 

e-3 


where  Hp  H2,  H3  and  are  the  thicknesses  of  the  four  layers. 
The  index  1  applies  to  the  surface  layer. 

Boundary  conditions  at  the  surface  (z  =  0): 

Tz^f  A,  ♦  B, -c,  C^.2(a,)  » 

X,»=o  A,  -  =0 

Boundary  conditions  at  the  first  interface  (z  =  H^): 

tx  •.  f  CT,  (jH-  2j^, -x)  e'^  +  J>,  (yt-  f-x  )  e- 

1-  33  -  Cl  2^,.x)  V  3),  [»-  2^.  e 

Tra;  A,^’* (^2yx,+x) 

Alt’*  -  Cl  C2vvti-x)e''  4  3>, 


w:  hX'-  c,  3>.  r 

’*-  Ca  --T>2  (2-4  ^■®  J 

U.;  A,e*  4-  B.eT’*  4-  C,  (><4-x)e^  -  2>,  (>i-x)e"^j 
Pw  [^/^2«-*  f  bzc'^i-  CzQu'x)^’^ -1>7. 

Boundary  conditions  at  the  second  interface  (z  =  Hj4-  H2): 

^z  -  Cl  (>»-2)A.i-y3%^  ^  Da 

T„:  C,ctyx*y)>y  V  Da C^V-A-y)^^  ' 

AsaV-  Bj.C’'  *  CjC2)A,*y)Ay  +  '3i»C^)A>-y)*'’’ 

1-..  [A,!”-  -  Cy  cs-^iAj-y)^’' J 

U;  AjaV  t  J,t-V  V  Cl  (Aly)ty  -  1li(A-y)»‘?'  = 

ki.  Mjt’’ t- 'bji'' y  C>fAfy3«.y -D>(A-y)*-’'^ 
Boundary  conditions  at  the  third  interface  (z  =  Hj  4-  H2  4-  H^): 

(Tz  :  A^«*4  Cj,  (>t-a^>.a)e*  4-  ‘5)3(m-2^3»-'*)^"*  ^ 

vC>4  0y^A  +  »)®‘*'  ^ 

t*.  •  C>  +  -  1>>  fA-x)<~'^S^ 

Lu.  4  6>^<r*  4^  Cz,^A4»)t*  -  ^ 

Boundary  conditions  at  the  bottom  (z=  Hj  4-  H2  4-  4-  H^): 


/'1. 3 


2.  Resolution  of  the  system  of  16  equations. 

In  the  equations  of  the  conditions  at  the  third  interface,  is  replaced 
by  its  value  taken  from  the  fixed  bottom  condition: 

We  write  the  conditions  at  the  third  interface  in  matrixform 

We  invert 

where  » 


nV^  - 


o 

-\ 
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o 

I 
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-\ 
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A  » 

- «, 


\  ° 

I  o  o 


0  1 

o  o  ® 

i  -V  ' 

o 

.  i 

1 

-Va. 

e 

f  1 

1 

1 

[2-4f.A4-V) 

-X 

4-  C 

1  )  } 

-lw 

Lw 
^  kw 

U(Z- 4^4  +  0  j 

'hO  L 


» 


I  .1  j  ^  a 

%-  vr  »,•  k  " 
'/-■.A/--' 

‘.V>.  A 
r. 

■■-[-‘•.■••.'A 
;y%v> 

S 

tv 

>  -  .V  •  1~i 

'* 

'■'V^V-V'V 
.  “  * 

■*-  •’  • 

>' A\  -■>. 

'-V-V-’S'l' 
•-•  -.'  %*  •, 

•;-<y  •:•<;••• 


%■'  %"  • 


I  r*:-.x’V Vvv* 


-  i  Iw  -v 

_  V  s  Uw  -  Lu.  ,0. 


WG  write  the  conditions  at  the  second  interface  in  matrixform 

CA2B.C-2  V)3^  M4  f/»3B3C3l>3)*^ 


-\ 


A[A']^2.) 


“V 

•  e 


j-[^yy)  -.(2,-Ay  ~[2yj^y)  -(^.2^,7) 


1 

0 


—I 


o 

~{a  -y) 

o 

-\ 


o 

-I 


o 

( 


1  0 

1  0 

Uw  0 

-kv(^~>^h3*7)  °  , 

kv.  0 

kwC-^^v)  ' 

I 

f  ' 

0 

\ 

f 

0 

0  12-^3 -y) 

1  0  -kv 

0  -Vw(}--^)r-i^ 

y) 

i. 

1 

“i 

1 

0 

A  ^  / 

-  f  n,i  •  e"^  + 

M 

32«- 


1 


T 


n>,.nAa 


O  V 

o  0 


O  v 

o  o 


o  o 


O  4- 

o  o 


O  o  o 

O  -I-  o 

o  o  o 

0  -t"  o 


^>241 


Tl>v^Al  ^92-T1a2 


-Xa 

o 

)<2 

c 

o 

ka  / 

o 

*  o 

o 

o  kj 

o 

O  0 

to 

o  O 

f*V3  o  to 

fc  kx  o  CV.-'^O 


o  to  VC^  o 

O  to  o 


V<,«  -  ^yt  ( kv - k 

V-i  e  — '1— 2Vcyy  +-<4  kw  .y  9  —  v»- 

k  j  a  kvv  -  kw 

—  ^)*-iK3  ”  ^ yt  k'Vyr  4  Syx y 3  Vy/  +  ku  —  ^yi ku 


■•■.  V  '.Va 

f"‘--  ■-'!>>.-/, 


'.-'/■v'-V. 


.V".\’v‘.'«.*  vv 


rj  - f*r-\ 


*'.*  > 


:'v'v1 


*  to  ‘  1.1 


-  •/  •.**  ' 
f  •  -  '  •  •  . 

•  •  .  •*,  to*.  •*..  •* 

■/  V 


nb.^1A  = 


0  ol 


Mi!>.  M?ju  -- 


«  O 


+  4^' 

I  C  0 

4-  4 

0  0 

i  +  + 


^  ^\Jk 


11 A 


^^f»(iz  =  ■>■  ■*’l  c  1^13  ^ 


hrb.  ^$>02  i 


r  IVi3 


‘.••V 


>  •  fc  ■  k* 


.  •  .'■•  .;  .'•  A 


-  /-  .*•  ^ 


M,,,, .'Mjij,  =  ■«  =<2 


(AiftiCit,)' 


^‘(''->'O(-'-h0 


We  write  the  conditions  at  the  first  interface  in  matrixform 

[A,rh,c.v,y,  n,  («.,!!>,  c^.J),y 


.-••y/'.-'-.-l 

.  .  '  A  .••  ■•> 

iJmJ 

MM 


(A.^.C.t:,y,  n;'.  M,  (A,fc,C,  Pt) 


^A:>:-:;: 


C<T)*)  «  -  — — — - -  T'  I 

4hi,.«-  4  li\.«.  .«.  t  lu-^  ^  ^ 

^  ~U  T-  ->  -2.(V--y)  -r  -^-b-y)  t-  ~%  -2.(.^ 


We  write  the  conditions  at  the  surface  in  matrixform: 


1  1  -  (4-2|>w>')  ^>1-2^*) 

^  -N  T 

c.  vO 

1  -(  2^1  1 

1  1  1  - 

C4-2j^.i‘) 

1.  -, 

2 

2^1 

4  4 

MI.Tij=  Uij  « 

V  4- 

In  .V  -il^-y)  ,1  -2.(z-y)  I,  -2lK;0,  ||  -h 

I  Uj)  +  0)1  +-  e.  4-  t/jjtj  a,  ■♦•  S'  •  S' 

^\1  ~^y  i\  "^-V-  1 1  -2*  i)  ->  1 1  _ 

4Ui2  «'  '  4  ^24®^  4  w^ie-  4  •*' 

,  V  r 


where  the  matrices  U..  are  all  (2x2)  square  matrices. 


We  write  the  terms  of  matrix  U 


Al.ll 


11 


u 


II 


ail 

O-XI 


and  the  sum  of  the  terms  of  all  other  matrices 


Cl, 

Ctt 

Cl, 

C-22 

We  notice  that  the  terms  a.,  contain  constants  and  linear  functions  of  the 

'  w 

variables  x,  y,  z  and  their  products. 

The  terms  c..  contain  all  negative  exponential  functions  of  the  variables 

*  J 

X,  y,  z  and  t.  For  high  values  of  the  negative  exponent  they  all  tend  to  zero 
The  matricial  equation  of  and  can  be  written  in  the  following  way: 


4  (cxii  4  c«)  Da  «  o 

The  solution  of  this  system  is: 
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V 


iA 

V 


4 


c.] 

C  11^ 


^  “  ^11  •  CX-Xl.  ~  C^li .  C»v»7  ^Ci  II  C  C  2 1 

-  [cKll  4  Ct|\  .Cii  4  C|,,Cl, 


It  - 


3.  Values  of  the  parametrs  . 


3.1.  Values  of  the  parameters  A^,  C^. 

At  the  bottom  of  the  first  layer,  parameters  Aj  and  Cj  are  factors  of  the 
positive  exponent  e^.  Thus  when  computing  stresses  in  the  first  layer  we 
have  to  input  at  least  two  positive  exponents  which  necessarly  will  lead 
to  overflowproblems .  Therefore  we  shall  express  next  modified  values 
of  the  parameters  A^  and  C^:  Aj.e  and  Cj.e  . 

Regarding  the  paramete"S  and  we  shall  replace  them  in  the  equations 
for  the  first  layer  in  function  of  Aj  and  Cj,  using  the  boundary  conditions 
at  the  surface: 


A,  +  b,  -  C,  (a- Ip..”) 

A,  -  1b,  C,.  2p,  =« 


The  resolution  of  this  system  leads  to 

A,  Ap.C, 

D,  ^  e, 


The  values  of  A^,  and  are  given  by 


^ p.) 


We  notice  that  the  matrices  Tj^.  and  T^^.  contain  nothing  but  zeros  in  their 
first  and  third  rows  and  that  the  matrices  T^^  and  T^^  contain  nothing  but 
zeros  in  their  second  and  fourth  rows,  so  that  the  values  of  Aj  and 
depend  only  on  the  matrices  ^2^  ^3i  ■ 


p 

p 


i T-jv.vr-"". 


'^.W'  *  ■.•*■ '.»  v»” 


■  '■  ■  w ■  » ■  u  "v  ■  *-'»  vriTB  V  ui"^'.  T  v^ '.■'  '.’* 


A1.14 


We  obtain  then 


e.  [  A,  o  c,  o)^=  - 


-^2y>x^  -liy-T)  ^ 

e  .  e,  V  Til.  t  ^  e 


-\^2V-?<^  _  _r22->')  ~  _x 
V  bA  "i-  ^ 


[t„. 

.y  '’‘-H^’-v)  y  -l[x~y)  -y  -y  -2lV-x)l  T 

We  notice  that  all  the  terms  contain  a  negative  exponential  function.  For 
high  values  of  the  variables,  the  numerators  tend  thus  to  zero  and  the 
denominator  to  the  constant  FKL. 


3.2.  Values  of  the  parameters  ^2^  ^2’  ^2’  ^2' 

For  the  reasons  explained  in  previous  paragraph  we  express  next  modified 
values  of  the  parameters;  A^e^,  B^e"^,  C^e^  ,  D^e”’^' 

The  values  of  the  parameters  are  given  by 


(A^  fc,  C,  ■* 


Ah 

V  J  .  [  bj,  3 


.  [v„ . 


r’t’-'V  V,  rZC-*'.  ..  Wn. 


We  split  this  relation  in 


y 
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(Aa  0  Cl  o)  c 
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w-c.-ivIn^-; 


3.3.  Values  of  the  parameters  A^,  C^,  D^: 

The  values  of  A^,  and  are  given  by 


- - r  ]  ^9101 


We  split  this  relation  to  obtain 


(A>,  0  Cj,  0)’’'=  ~ 


4(4-)-^) 


(o  t  ^ 


^  "  AMI 

3.4.  Values  of  the  parameters  B^,  D^: 
Following  the  same  procedure,  we  obtain 


r  ^ 


4.  Determination  of  the  stresses  and  the  disolacements 


The  program  calculates  the  stresses  and  the  displacements  using  following 
relations: 


a 


(Tz®^  I  3QC»*^r).  y 

“  V  -  Vm-X  -1 

_  _  vjL  1  c>^o .  j.  (‘^•0 
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Stresses  and  displacements  are  calculated  in  a  system  of  cylindrical 
coordinates.  The  stresses  due  to  several  loads  are  to  be  added  together. 
Therefore  we  must  express  them  in  cartesian  coordinates  using  following 
relations: 
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wherein  the  signification  of  the  angle  «  is  illustrated  below. 
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j  T»  Cv^r^ 

-  [A,t,’‘)(-^-2iwl»)  e  ^  \  (C,fc^)[Z)A,‘V>K‘li--4^,vt^4<.^  Jduu^ 


.^v 

1)0..  [vvifrv^ViAVe  cA.H*^r 


o) 


1^\  - 2—  .  K;v.>dL4  (r#.) 

•'o 


k 

lihl  ^  0,  Y  C ,  t 

-  T>,  J 

2f^-vn.  f  7.M.j,^»>i.--hj...  . 

E,  1  *"  ^  J. 

J”*  X_W^k)  J^(A,i')C^’‘'’“''l  (C,  i")  \.2-4t.,-'^'')i' 

I  )  Wn 


-n\  ^f'.x„^^^vr’-W'>^^''^ 


('('H 
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APPENDIX  2 


Algebraical  Analysis  of  a  four-layer  isotropic  System 
with  fixed  bottom,  full  slip  condition  at  the  first 
and  second  interface,  full  friction  at  the  third 
interface. 


APPENDIX  2 


Algebraical  Analysis  of  a  four-layer  isotropic  System  with  fixed 
bottom  and  full  slip  condition  at  the  first  and  second  interfaces 


1.  Boundary  conditions. 


We  write 


/A  i  -  A;m^  1>'\  = 

Ci=c?iv^ 

D  -  ^‘iWi 

c_ 

1  -•  - - - 

ti  0^ 

X  c  H , 

ym.  Wl  Vl  1+  H  x) 

2=  wCHii"  vi2'f^O 

V  -  Ml  ^3 

u  =■  -  a  \-+  z 

where  Hj,  H2,  and  are  the  thicknesses  of  the  four  layers. 
The  index  1  applies  to  the  first  layer. 

Boundary  conditions  at  the  surface  (z  =  0): 

(T*.).  A,tb, -C,  (^-V)  •  = 

Xri  =  <,  A,- »,  +  ^,.2^.,  + 

Boundary  conditions  at  the  first  interface  (z  =  H,): 

V  bj C,  4-3)i 

Tfz,so;  b,t“%  C. 

r„,=o..  -  Ne'  H  ^  3),  -o 

-C2(^Z-Aw,-;^y-  J 


,  ^.  J*:  ■ '.  '^ 


A2.2 


Boundary  conditions  at  the  third  interface  (z  =  Hj^  +  H2  +  H^) 

(Ta :  Ai*  f  hijC'-  Cj  (;>»-2^^-z3  J"  V  3^  (>>-  2)^5  r)  e 

"bjyC^  -C-A  ^  1>A 

■^fZi  ^3e*  f  C2,  ^3C^>'2.3e  = 

W  :  A^e?-  l>3ft:^-  C},(z-Ay^'  ^)^  -  5>iC2.-4)^5t*^)e^  “-- 

I J/  4 «?  -  c*'-  ca  (2-  -<  r  A  -r)'-^  -  ’’•<  {-2  -  ^  r-*  ^  'J 

U.  ;  >!,/<-  0,C*  A  Cs  (lA^^)*-*-  (>-^K'‘'= 

l  AAe.^ QC'+O*-^' t”' J 

Boundary  conditions  at  the  bottom  (z  =  Hj  +  +  H^): 

Cy,  c  o 


Boundary  conditions  at  the  second  interface  (z  =  Hj  +  H^): 

-Jtt — 1 

C-v  •-:-••  :.--:j 

ir^:  A2«,^v  B)7^^_Ca(>i-2^,-v)«.>'+ 

3)3(-i-v,+v3«:^ 

rrai=o.'  A  +  Cu  +-3i^2f<i'y  )«- 

A)*y  -  &»,«.*'^  3  Cs,Ca)-54y)«,y->- 

VY-  Az*? -  8>j^e.''^_  y)®' 

W  ■  *  •  V.  *  .  ►•  fc 

V‘^-  C3(2-A^3'vK- 

■-  >■  •  ■  -1 


.,-  -ri 


m. 


•*- <■ 

i  i: 


i  ' 


.VA'.-..-..-^ 


2,  Resolution  of  the  system  of  16  equations. 

In  the  equations  of  the  conditions  at  the  third  interface,  is 
replaced  by  its  value  taken  from  the  fixed  bottom  condition: 

We  write  the  conditions  at  the  third  interface  in  matrixform 

Me  ibA  y>Ay 

We  invert  \L 
b 


mV.  ibA  ^Ay 


where 


n*5  =  -  — j - r  ^ 


-0-2) 


I 

-U  > 

L  LCl-4y.uh) 

\l 


-7. 

Y-  e. 


C3-h-*) 

--L[^y 


_  .  M\5,  4^  t  (Ba! 


^  ^5tG2  -^^5241 

1-  C  ^Ay 


K' 

r 


J.»  ■!.'<  j  ■■  »■.'  ■■'■  ■  ir» 


We  consider  nov/  the  boundary  conditions  at  the  surface. 

Adding  and  substracting  the  surface  conditions,  we  obtain 

2  A)  *»• 

“lb,  =  X  -t-  C,  -  Til 

Adding  and  substracting  then  the  first  tv.'o  conditions  at  the  first 
interface,  \ie  obtain 

2Ay  <■  l-'O 

2i.t"  -C,..,’'  ,  t), 

where 

We  replace  A^  and  by  their  values  in  function  of  and  and 
solve  the  system 

C  e.’'-  -  (x+2x)  eT^ 

3)  e*  -v  (x-2x)  ^  _  e 

We  combine  the  T  -conditions  with  the  w-condition  so  that 
rz 

and  with  T,  =  F  _ 

C.e’*  ^  :=  F^  V 

We  replace  and  Dj  by  their  values  in  function  of  A^,  R2»  ^2  and  Dp 


■*-^x 

_  2(x-»-x)  «r^  4  2  1 

N  -2.H  -/t-A 

je  _  e  -  X 

-*•  N,'*  A 


Together  with  the  T  condition,  we  obtain  then  the  system 
^2^^  -  r  C2  ^  «.  -o 

T  bit"’'  -  Cx  ^  A-  D7(>1-2yiiTX-R5)t  r 

where 

A>r  A-aZ^" 

2  -  aC^t’<)e’' 

V,  =  2[Aflx^)Z^-  -i 

R. 

We  notice  that  for  m  =  *0 

Ra=o  7i^-^  R:>=  -‘F^ 

We  solve  the  system  for  A2  and 

2Ai*  -  lie’' 4-  Ca  (^->4^-^’' -*-^0 

?>%2  •  -V  Cl  -  Di 

%y 

Adding  and  substracting  now  the  first  tv/o  conditions  at  the  second 
interface,  we  obtain 

lA^e^-  Cl  A  )m -2y  ^  •»■  5i  e  ^  =  ^^^3 

iW  -  C2e>  +  Di(^-/^i-»2y) 

where 

[/\b]  =  ^  5,  (^-2^3^  yK^ 

We  replace  and  B2  by  their  values  in  function  of  C2  and  02 
and  solve  the  system. 


A2.8 


■  V,  V  W.  V  V. 


We  notice  that  for  m  =  #0 

We  write  the  system  in  matrix  form 

-  [A-2^>-y  vVSj)  o 

O  '^'1' 


)  o 


o 


-I 


'l)-h-y 


- 

*} 

1'"' 

- 

t. 

-h 


or 


y 

y.  ma,  ^  ’»» i'^^  *) 

We  replace  the  matrix  (A^  in  function  of  and 

c  ^  ■*) 


We  notice  tliat  ~ 


^41- ^'5262  "  ° 


^2-^161  =  0 


^42 -^’5162  "  ° 


*  %*  ST 

I*  ^ 


^  *  •  *  -  -  * 


m. 

.  ^.  " 


N. 


V  _"v 


•  •  •  c  •  •> 


,  ••. 


«  •  f  •  V."  • 


v*. 


f  • 


V7  IfS'Vv' 


VJe  call 


^42‘^^!>261  "  *^3 


^-12 -^^'0262  " 
and  write 

[e+'-'>  w,  V  ,  .-’'v*J  Lb*  M 

=  >^(-'-h»)(  %  “3”^ 

that  we  transform  into 

We  write 


-(NV)  ,  , 


Al, 

ait 

ft  J> 

a  tt 

IVi  = 


W> 


where  all  the  terms  a.,  converge  v/hen  ni  = 

*  \0 

We  devellop  the  matrix  cqiiation 


(o.„  V  1>„)  l*er’'  t  («.v  *  t,.  )  P*«  ’’  •  -'/-'-(■>)  ^ 

(,„  ^  u.)  hiC  ^  [«.*  ♦  V»  )  '  “ 


We  solve  the  system 


A2.10 


*-l2  1*1 

*.  tVi.  4-  Wl 


where 


V^=  ( cx,x  »- lo.OC®'*-J-^  '•^O  l>ir)(‘ttj.i  4^  Wi} 

=•  [^CX„4  Uli)  CVj.2  -  (^CK,X  ♦  ^  hu  h  12  ~  h)x  hz\  4-»-W  i>12 

The  term  ^11-^22  ~  *^12‘*^21  linear  functions  of  the  variables 

and  has  to  be  devel loped  in  close  form 

I  0  0  I 


'on 

1  '  ' 

0 

111 

9k 

biz 

1  “  - 

i  0 

1 — >1 

\>I7 

=  Li,4-X 

U,  (A. 

a  -L, 

Wz 

= 

-aCU.-Lz) 

»■  L, 

C2-k>'y) 

bn- 

bzz  ~  ^11- 

bjt  s  L, 

I-lv 

.  X  (l.-Li)  i-  1- 

4-  4 

- 

-LO  V  u 

=  L. 

Lx  - 

-c&>) 

The  linear  functions  of  the  variables  have  disappeared  so  that  for  m  ^ 
the  numerators  of  B^.e"^  and  D^.e’’^  tend  both  to  zero,  because  of  the 
factor  Q^,  and  the  denominator  tends  to  a  constant: 

)iw^  <5,  L|  Li  k-))  I-.L2 

and  finally 

IV  -X  n.  .-y  rf*'/)  -K  -I  r.  -y 


3.  Values  of  the  parameters  ,  D^, 


We  express  the  values  of  the  parameters  A^ ,  D|  in  the  same  way  as 
explained  in  appendix  1. 

3.1.  Values  of  tl>e  parameters  A^,  C3, 

The  values  of  tt\e  parameters  A^,  B^,  and  arc  obtained  from  tiu' 
relation 

Ci,  ..  j  1 


The  matrices  f1t326l  ^'5262  nothing  but  zeros  in  their 

first  and  third  rows,  so  that  we  can  v/rite 


(Aj,  o  c>> 


DO 


^  -[V-y)  ^ 

- - 1.  .4  ^  >  '5\i  I 

♦  (^e’’ 


The  matrices  and  ^5^52  contain  nothing  but  zeros  in  their 

second  and  fourth  rov/s,  so  that  we  can  write 


(° 


.  -V  .  tm-t  _  _  J _ 4  otitiT 


r 


3.2.  Values  of  the  parameters  A^,  C^,  D^. 


We  have 

C^t.  -  - - 

SO  that 


where 


We  have  also  that 


2Ai<t'^-  Ca.  V  s 


so  that 

A,.,7=  J_|YA3]  t  C,(vAv..-27)t’'  - 

and  finally  we  have  that 

"XYfi  s  —  Cl  ^^3  ~  ~ 

Hi 


SO  that 


A?.  13 


3.3.  Values  of  the  parameters  Aj 

and  fj. 

We  have  that 

«  [All  y  ux-.K>'.  .1 

_  2x)C*  \  €^'‘ 

,  e  - - —  ■  ■  ■ 

V. 

V. 

where 

We  have  also  that 

so  that 

ij  L'^O  -  c,[^-At.,->-)‘'’‘-  ■^'*  ) 

The  values  of  Bj  and  Dj  are  obtained  from  the  surface  conditions 
A,  >  -  C,  4  5),  (-1'  ■2^.)  =  ■< 

A  ,  -  lb,  r  C,  2y.,  i,  2^,  =  o 

so  ttiat 

3,  =  -/I-  2  A,  >  c,  C^-A  )^.) 

s  2(A,t’')e"’‘  V  (C,  e’')  ('l-A^'.)  fc 


4.  Relations  for  the  stresses  and  the  displacements 


The  relations  for  the  stresses  and  displacements  are  completely  the  same 
as  those  deve'lopped  in  appendix  1,  by  replacing  the  parameters 
by  their  adequate  values. 

Nevertheless,  there  is  a  problem  in  the  computation  of  the  vertical 
displacement:  its  value  at  the  origin  (m  =  0)  is  undeterminated. 

The  relation  for  the  vertical  deflection  at  the  surface  is  given  by 

^  [  A.-  B.  -C.  -  ■>, 

T  j  L 

which,  to  avoid  convirgency  problems,  is  transformed  into 

The  numerators  and  denominator  of  A^  and  Cj  are  both  zero  for  m  =  0. 

To  eliminate  the  indetermination  we  should  develop  Aj  and  Cj  in  a 
Taylor  series.  Altough  this  is  theoretically  possible,  the  required 
computation  is  very  long  and  the  risks  of  introducing  errors,  in  doing 
so,  are  enormous.  Fortunately,  we  dispose  over  the  fixed  bottom  con¬ 
dition,  which,  for  m  =  0,  transforms  into 

The  w-conditions  at  the  other  interfaces  transform  into 

A-,- hr  -Da(2-A|..)=  W  -Ci 

A,-  lb,  -  C.  -B?  (2-A^<}-  F 

so  that,  for  m  =  0, 

A.  -  t),  -C,  -D, 

and  thus 

J^A,  -  =  'T 


so  that  the  problem  is  solved  without  any  difficulty. 


APPENDIX  3 


Algebraical  Analysis  of  a  four-layer  anisotropic  System 
with  fixed  bottom  and  partial  or  full  friction  interface 

conditions. 


1.  Boundary  conditions. 


Because  of  the  presence  of  4  different  exponential  functions  in  the 
expressions  for  stresses  and  displacements,  the  algebraical  analysis 
is  more  laborious  than  in  the  isotropic  case  (appendices  1  and  2)  and 
has  to  be  developped  in  a  more  detailed  way. 

To  reduce  the  number  of  exponentials  we  write 


x  X  f  •  M.)  . 

VV(5<  H  ,+  Vlxh  ■  HO 

C'^»v^  XiSiQ^i+  y'O  ^ 
iw*  n\  Si  (yii  +  ® 


Further  we  also  write 


X  - 

y  r  (  VI  1  4 

Z  =  VV,  ^3") 

W.  (W  ,  K  Ht  4-W3  ^ 


where  Hp  and  are  the  thicknesses  of  the  four  layers. 

The  index  1  applies  to  the  first  layer. 


Boundary  conditions  at  the  surface  (z 

(Te  =  ^  •. 


0): 


-X 


_  B,  y.  C,e  ’  -  i-.li,  =  o 


■  W  ^  --n  L'>«  fV  ^  \  -V-  ■C’’  H. 


Boundary  conditions  at  the  first  interface  (z  =  Hj^) 

Tz  :  /i .  b.  C,  J),  e^' A 1  E)a  +  Cz  ^  D, 

Xa-.  fte’'  +^,C.-i.]>.e''^=  b^  -  s.  P, 

vv:  (^4yA,)  A ,  -  C-A +p-i')  ■*■  *^iC'^'+^0  ^  ~ 

F,  A  jc. —  t,-^-*-F' ' )  ^2  +■  J I  ^  -  5i 

U.(1^.4F.)A«  f  B.  <-■’'■»-  L -14- f-.)C,  4 

Ti  j^[V)afF‘v') 

Boundary  conditions  at  the  second  interface  (z  =  -f  H2): 

Oz:  C,  V  A3e-^"->\E>3  *  1>3 

Tfx:  Aj-  SiCx  -  b3  ''1  x»  J) 

W'.  ('1+f*^)Ai  -  (’-(->^2')  i7(>>*^/^2')Cj  +  ^■' 

(44^-4)  Bi  4  ^0 

U:  (h7tjA2)/,  V 

Vfi[_(4'!.-»F^) Bi  +■  J 

Boundary  conditions  at  the  third  interface  (z=  Hj  ■^  H2  +  H^): 

L,  J(-l4F2^)A>;e^^'*2  (>i4.^a)B/,  >4  ^-S^ihA^-y^)  C> 

L7  ^  *-(^  +  h^)h4 

Boundary  condition  at  the  bottom  (z  =  -4  H2  +  H^  +  H^); 


tv-.  (><4|aa)Az  -  (-4-4^a)  B4  *  J/,(h^A*f/yA)^A  - -5^4  (t^A  4^|*A)iU  «- 


-.J4|V'  = 


2.  Expression  of  the  boundary  conditions  in  matrixform. 

2.1.  At  the  third  interface. 

In  the  equations  at  the  third  interface,  A^,  or  C^,  is  replaced  by  its 
value  obtained  from  the  fixed  bottom  condition. 

We  write  the  conditions  at  the  third  interface  in  matrixform 


We  invert  M. 


03)"^=  Mi'. 


•^1  A 
I5  - - - 


■f  s.S  1 


-2(_V-al 

I 

-1  +  a- 


f  S4<1 


1  + - < 


U.y4>>*4Vi,)L-^  ^  ^  J 


214  {y-») 
-Zs^lhfL) 


We  write  Cj  and  in  function  of  and  D^. 


^  ^  Y^(.b^)  hA  V  ^  HTCa;)  .  V>fcf>2)  Pa  J.’^-y) 

j) 

^  :Lr^1(^^).  rAiiu).  \b4  -t-  ^  M5(3«)- 

^3»  =  - - - 

_  Mc/i>i)- bi?  +  ^ 

U  ^  =  — - — — -  c. 


!)>>  = 


where  Mg(i,j)  are  the  constants  in 


We  write 


T>i^  = 


m-j60  n  ccm) 

M»(jO  nt6^) 


vv.v.!l*.:Sj 


.v.->.r4 


• '  •-*  •< 


•  /-v-.v. .. 

'  y  ••  W  •  .  •«  fcN 

•:<-‘.>.v.- 
-  ••  ■ 


^  ’«*•  .  ■  .*“  - 
w*  '  •  •  ‘ 

^  r  •  La 


so  that 


P/.,e 


Ja{»-v) 


j,(a-y) 


2.2.  At  the  second  interface. 

We  write  the  conditions  at  the  second  interface  in  matrixform. 

Mj  (.^i  t^h  ^3  ^3  P3) 


l/e  invert  M. 


.'.VJ 


l5t(vl-h*) 


^„J)-')  ...Jv*) 

-'  St 


-(a-v) 


-4(t-y) 


K(^+|»'s]c  ^ 

Wi(m+K')  ki 


-v<,r>K4|^i) 

Vf^  (4f^j) 


We  write  B2,  C^  and  D2  in  function  of  A^,  B3,  C3  and  D^. 


5tXi  I'X  -v»i) 


5Ji  ( -4-  v»  t') 


^Ji  (A-v^a) 

2.  Vi^C^O- ^13(30- V)a /  iA)  1>3 


Wl  (A-V\x) 

-s:  V  ^r^>(A\).'nAM  ^3  /Wv-^') 

J^3»(a-v\  i") 

Pa- v< 7^ 


A3. 6 


where  and  are  the  constants  in  and  M^. 

We  write 


so  that 


Aa  ' 

hz 

A 

c. 

D.1 

2.3. 

At  the  first 

Q.t*"'''* 


«>> 

C> 

Dj 


We  write  the  conditions  at  the  first  interface  in  matrixform. 
[A.lb.C.  D.)'r  =  Via 

We  invert 

(  A,  b.C.  •&,)'’'=  (A,  ftiCi  ^2)^ 


m:’-- 


Ha* 


-y,  ( v>,+  y^«') 

5^  ~  r, 

A 

jC 

-Ti  (v».>|^') 

-J'xCV'v+V'')^ 

( 

r  N 

-» 

1  1 

\ 

-St 

We  write  Ap  Bp  Cj  and  in  function  of  ^2’  ^2  ^2 


A3./ 


A, 


C,  = 


SS) 

aj-, 

T),  -  - - —  .  e. 

ii,  ( A-y> <) 

Zs,  i^-y^  0 

Q  ^  ^  ^H^[30.  V>^6^)  £>? 

as",  ^A~  v> 

5  H./-3  0  -»  s:vl-,c^iyM^o^)P2 

ii",  (a~y\  .') 

T)  „  ^  ^H,(4>)-n£0^)S7 

£  M.(4i).M,(n)«:’'''''h  »  s-yi^i^‘).n>0^)^2 

4  _ _ _ — - - - - — ^ — — - —  -e 

where  Mj(i,j)  and  M2(i,j)  are  the  constants  in  and  M^. 

We  write 

Rj-i=  M. 

SO  that 


T>,= 


p  rh-*) 

K\\t. 

^lit 

Rne^ 

ftv.i''-’' 

(!«.-'''".V'’ 

0  -;,tv-*)  /.X 

t 

^44 

•.••‘.vv.'/ 

•\  •*•  -'.  •‘‘.  •’ 

.*  ^  \*  *.• 


‘  ^  0*0 


,v.'*\-V-V 

f 


y:-.-r;<v:- 


!•  'Sv-'N' 

-''V*  .vV," 


^Resolution  of  the  system  of  boundary  conditions. 
We  write  the  conditions  at  the  surface  in  matrixform 


-1 


-r.y 

e  4 


-■A  s.«.  - 


•  (A,  C. 

We  have  then  following  system  of  matrix  equations 

hr .  M.ib.c. 

(A.lb.c.  _  HR 

2r,f>t-v>0 

-  M(^  (A^B3C^ 

Jri  i") 

- MP  i  6/ 


^^3  (>»-»' O 


SO  that 


I 

o 


- - -  .  Hr .  nY- 


Hr.hR  = 


a, 


ail 


V>i> 

\>ai 

b/ii 


tXix  CK.\y  cx »/, 

(XXI 

ba 

bii 

V>3t 
biiz 


Cl,  0,1. 

C2, 


rir.MR.HQ.HP-c 


cir 


s  ?  Ji  Tl  *•>  • 


Ci-l 


Cji-  Cii  -  C  II  .  Cli 


Da  =  -  rjr>, 


_ C^t _ 

til  ■  Cj j .  C.I2 .  Ca  I 


3.1.  Determination  of  ^^^022  "  ^21*^12 

The  expression  of  the  denominator  ^^^.022  "  established 

in  complete  closeform  al tough  regarding  the  exponentials. 


Cm  = 

0.1,  -bii 

^  o-ix-bj.,  + 

0-i>-  bj, 

CX|^  .  b.<,i 

Cix  = 

cxn  b,2. 

V  oiii.  bxx  + 

onx-bji. 

¥  C\i4.  b.4x 

Cx,  = 

Oil-  b(, 

4-  O.XX.  ba,  -f 

¥  CXj^bAi 

Cix  ^ 

KAi\  bn 

♦"  aix.  bii  ♦ 

CTij.bj? 

V  A-3A  ■)>A 

Cii  Cai  « 

{cx^^  eL^^-  ati.txii)  (\>ii  bix  -  W 

f  (wn  ai3.  cii»y  (bn  bn  -  bu  b#i) 

f  (cxii-  -  oiiA-t^'O  • 

T  (cxa.  ax^-  A,^.cx,2).(b.,  bn' W  b,0 
V  (cMx.aj/,  -  ajA- 

¥  (  oii^.cilA  -  PDA  b>i.  bA\) 

] 

] 

] 

J 


a.,a« 


_/.r 

-Mf 


Rai 

.  Xi> 


CX|A=  L. 


I^3A\ 


[-a1^23  f 

+  Ra2|_a<^»2.ft:V^'’'  +  •*-  f^'J')R^3  j 

^  ^/a  j  ♦-  e»‘‘'')^i^'  V  21-, 


S^yt  y 


*  Sx»  [  2*u  ^  C-*')  3 


+  R  32  ]_  C^-.^')  ^tN  «- 
4  R^2 


^  -S-/*  1 

^  4  '2r,  e  J 


-jjty-*)  -^.x  x’ 
tM>!»  -  «-  e  •  ■ 


-2J)* 


r  -y  -J-Iif  4.  U+JO^JY  t"  '  4f'<-^'^^YvJ 

4-  IZ,  W  K>A  ^  .X  ~J,v-\ 

U.3  ,  -y./.yl 

,,,,  .-.A—  J 


-IV-y)  J,«  ^.<  4,^ 

0\)(.  OXll  -  0'U  <^y4  = 

Ayi> 

OX„.  c^l^»  -  ^  ^  . 

t^-y^  J,y  A»i 

ox,,.  ^ 

^-x,l7-y)^rxy  ^y 

Wn.  C«Z  J  _  dcx^.  CA2T  - 

^  ^  Dxii.  DX51  5  «.’'.  Aa/i 

C^ia-  ^2A-  CIlA-  ^22 

-5,t7-’^').  J'’'.  e.^  A>2, 

w»-^.  0»  24  -  ®aJ4  -  *- 


/C\  O-  V  h-L  *■ 


a,  ^ 

>n  --  [l^2v  ^a  ^ 

b>,= 

b„  =  <l>,.Pa  •■ 

K,  =  ‘  ‘ 


^  ^72-^2\ 

4.  (PlX-l’l'i 

4  ^n-ht 
V  ^>2-^1^ 

4 

4 


4  (&iX^32 

♦  (&2j.p3. 

4  <Jl2}.P)2 
♦  ^33- ^a» 


^AJ-^32 


)  rt  D  ^3  (*■'*') 

4  lS)A 

4 

4 

<  «♦*  ’I 

i  Hm  f'lj]**'’’ 


V  ^,  .(V-1 

.'bli  -  \)n.b2t  =  •• 

4  'll!!  -  4l>  .l{j-v) 

4  -  C>\l  L^^''  .Jj.y")  .X>(f-V) 

v’[4,v(>i‘.-  «lt>]l’’»'  P«- 

,  [u..,®..- 0.4 «»]!:?«.  P4.-?«p-j ) 

„  ,  -2f2-y) 

V„  4  t»...  P'^  -  '»•  '■’0  '•  ' 

‘1  ">  o  "\  -^»C* 

=  f  ?W-  ^^2  *' 


\v.v.'.'> 

'"/•^.V.v.v 


%  ,%  %  .'X 

■■v-y-’/-  ►'• 
'■^■'.y-’.'•.^ 


'.■'.v-.-f.'-* 

k  •  »  -  fc  \W  •.• 


I  .  t  ,  «  .  M 


'.vl-.'.s'.y 


4^(£)a. -  ^)r»  <D3i3- ^U)>j 

^  ^)/4  -  <D\\i  *C)33  ^  ^32  "  ®)2  ^03^J)»'^^j 

.  V  V  U 

^w-  ^  om*  bA\  «  ^  ^  ^ 

I.V^xQAi-<OR  OA.]  >>a  [v^«.<0aa-  <£).a-C4.]  |>.a 

+  ['0)i-<JAb-  ^>3  *  [(^IbQl^t  -  '0l>  (tA3}  y>»» 

+  y<0l^  ®AA-  >3A  3-  [{^ia(S)A2  '  Q>il- 

,  X  X  In->')  (*-v)  e-*>l*“>') 

Wt  Wl  -  ^2*'  re® 

j  l-ShvOn-^n  tj']  )).^  ^  [«,„.<l!.A- 4)J*  »)'3  V'* 

f [vSivOjA - «»■•  ®'0 >»*  *  l''i>*'*”-  i 

V  V  N.  >. 

V?!.  \)A2-  >)JX  OAI  -  «- 

hit.,  w  ‘’“H*"'' 

'  ^HiiHai-  <>ii  fl4>3'”*  *  I’''**'”'''  ■  '” 

*  "aA  -  »1A  *'*’3  1^'  ®“'  y 

V  V  t,  k  ,--'»-'h<’-’'h‘><'-"), 

\)3)  DA2  -  P)2  ^Ai  = 

®»j- *  i^hA" 

A  \  «AV«AA  -  <»,A.»AAyV>A  <  t®>A.®.A  - 
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bn  ^■»2  -  = 

)>\\-  "bAx  ~  "bii  -  b/|i  = 
tx).  Wi- 

\'i\  \)A2  -  bxT.  bA>  • 
bii-  )>Mr  - 


j2-v) 

2>,3 


The  value  of  the  denominator  is  then  finally 

V  t.y  i»tv-^>  „C*-y)  ^^»C*'y) 

C|,.  C7X-  C,,.C2,  =  .  c. 


Ao.  ^ 


-tv-y) 

-2('#--y)  -iJiCv-’') 


+  Ai2).  ^>3  •  ^ 

»  A, A.  ^\A 

■V  Ai3  •  ^23 

1-  Aia  • 


-tv-y)  -SiCy-y)  ‘\ 

^  A3A.  ^  J 


The  term  contains  the  constant 


1Q12A  <^A1-  ^n4AAj  [^.-^22-  ‘^A7.^2.] 


We  v.rite 


C|,.  C  22  -  Cil.  C2,  =  «.  ■ 


X  j,y  (v-’^  C 

•  Q-.C-  .W  -c-  •  ^ — 


^2 


P3A5 


3.2.  Detarmi nation  of  the  parameters  and  C^. 
We  write 


\ 

Ci^,  = 

L^h- 

-lx  -J.< 
«.  •  e 

-  V  ft-  • 

\ 

CX27  = 

-1a  -r>x 
e  .  e 

-Vf.  4 

1 

DVlo,  =. 

^TX  -JiX 
t  - «. 

~  R 1"^  C-  '^*  4  J»  R  J'i  t 

lft.v 

-In  -J  .X 
fc  .  e 

j,Rj,v  eH 

b'n  » 

[4,. 

D 

•  <-  •  < 

-v).  0  -(a-y)‘ 

b\i  = 

^.7 

^  •  R)!  * 

P>7 

+•  CKv  .  Pi  2 

w,, . 

^  Pi» 

♦  (Dm  •  Pii  J 

Va.  , 

1 

Ja 

^  ^)22-R7> 

»-  ^2^  Pi? 

*  <Ch  fi,  J 

^\\  r 

L. 

[Om 

•  P., 

<■  (J>2  R» 

h^ 

■*  Cij^  .R*,  1 

[Q, 

,.Pu 

V  <0}2-  Pl2 

+■ 

J 

V.,  . 

r 

»-  Om 

+  P>. 

i\r  ^ 

4-  (J;,i  Ri 

+  (^44.  Pi? 

4  J 

We  have  that 


lb/,  =  %S',S'')  , 

D/<  -  ~  5  r.  r  3 


Ci2 _ 

c»».  CX2  _  CiJ.-  C-  li 
c  21 

C),  CX2-  Cil.  C2I 


!•> 

k-s 


* 


w^^^9W9m 
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4.  Values  of  the  parameters  A^ ,  D^ . 

4.1.  Values  of  the  parameters  A^.  B^.  C^,  0^. 

The  values  of  the  parameters  are  obtained  from  the  matrix  equation  in  §  2.1, 


A,  1 

ih 

V>.7 

bJ 

A 

ft,  el*-) 

9 

*77 1- 

Cj 

hi 

ftj 

One 


obtains  immediately  the  values  of  the  parameters  A^  and 


As  = 
C^c 


[>»- 

A 


+  ^17 

[Pv.  B.(  <- 


The  determination  of  the  values  of  the  parameters  and  need 
some  more  computation  to  insure  converqency. 


^3' 


•  \  ^  ^  ^  ^  ^  J 

__  \  P^,.  B/,  4  P47.  D/,] 
1  -* 


Those  relations  contain  positive  exponents  which  much  disappear  to 
avoid  overfh^w  problems. 

eC’-y). 

\z 


^27- 


T)3-- 


r  -  ■'*'  ^  *■/ 

K‘-WAy- 

•  .■'  ."•jr',’' 

»  .  ■  .  '  .^  . \ 

-V 

*  .*'  .*■'  •“', 
v' 

*•"  %“  o 

-  *•  i^,** .  *  -  • 

- 

t 

f.. 

XX' 


.  C\2 


VVil 
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?J,.  6*  -  oi  = 

k  «->’(>-'■>  ^Pj,.  b'oi  -?al>’>] 

^  kl„  .  ?„.  Vj.] 

+  A*  I  P2I 


-922\*»  p..  ^ 

.e^’-'’f22[(i..  P''  ‘ 

=  e ■  "Pfew. 

P2V  ^'12  "  P2I  ^21  - 

^-I’-V)  P,,  ^2,.  SWi'”'''*-”'*-’'’  ♦  O2V  ^22'“'”''’  *  H 
-  V22  t«)2..  P«  "  ®  »  *’»'  *  ®'* 

.  g-l’-V).  ?B22 


v>j 


'.V.^- V.%-J 

■,  •-■.'-■.••'.VJ 


;-:.-»r:-s 


k ».  -  ,  %  » 


'  W ''  -  •  '1  ^ 

.y>2-\y.V 

-V»'.  -*  •‘.  . 


;.  ^  iK  i 

••’/avL’.-!'-', 

•  •  ?  ^  • .  - .  • 


.•w>v/.v.-3 
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?at.  = 


_  ?7a  V„  -"  ^^5-?3)  ^ 

'?2\  '^\l  -  ^-*1  ~ 

^-(l-y)  Pi,  |^Q„.V„e^’'''V”'”''  *  1^4>  '’32e’’‘*’'’'4 (»*;,■*’«] 
-4.-^’‘'*.P.l\^0«.P.,  4  «4j,P„  4  44,.P4j 

,  e  <’■''’  ^6^^. 

The  positive  exponent  e^^  can  now  be  eliminated 

b3=  >^S,  J-a(»'V»0{^>l''^l)  • 

r -2(v-')  -t.(y-«)^'^_  'Pft<4  4  <r'''^''V22.  p«>« 


4  ^-W-’'’  -J->b-')  -p E,,J  4  «'24 .  T 0^  ]• 

and  the  numerator  contains  again  only  negative  exponents. 

Tv,,  4  C  ■4WV.2 


C41 


4  ,  e'‘'*-’‘''c4'24.b'4,^ 


-ft,4  U"''''’‘'e-”‘’'’‘’r,'2..b'n  4e-”‘''-'’cW>'v, 

^  -2iv-4V>4ir>'V„  [P4,.b',2  >>’.,] 

4  ^lr’'V”’ll-’'''  ^P4I. b'w-  Pm  Vj.]  4£^’‘''V24 [0,,. b«-  3 


^  '.•  *>  • 


/«* 


yy-y-y- 
.  -■.  ■■  •■•  •'. 

•■■.••.-  •  /  > 

av^:: 


-,-W  / 
\'  .•  -*  .• 

.*• 


•  •  ■•''•"  N.  *. 

•  •  *  •  • 


>  /•  .  '  '  V- 


’*  **  -,*  V  s*  ** 
•  •  ■  .**  .• 
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V>\a  -  =• 

-Pm[u»Pi.  ] 

_  t'*'''  ’’!  ? Ill  ^tSi.  i’ll  *  ®ii  Pii  *  'Oil  J 

,  Jilvv)  P6/.1 

?A).b31-  Va1>5I  -  _|4.y)  ./i(»-V^1 

_  j- p„i’l’-’')  »  5)„  P„  +»ii  *•  -i 


y^J,^ 


?A,  V*i- 


_^..,\aw)  p„^on,  P„ 


-la-v)-»>|»->’'l'\ 
.  iiAl  *-  ^A>  P32t  ^  J 


~ir>  li'V) 
=r  e 


Vb  Alt 


and  the  positive  exponent  e^3^^  can  be  eliminated 

p2,=  4  s.  ri  • 

^  ?eiAi  1-  "'lA- 


D-^a 


4.2.  Values  of  the  parameters  A^,  62 >  C2*  O2. 


The  values  of  the  parameters  ^2  ^nd  C2  are  immediately  obtained  from 
the  matrix  equation  in  §  2.2 


is, 


A,  t ^  V  iSJiA-CiJ 


c,  .  -  r®j,.  »  cS)5v6j 

2j,  h-Hi')  ^ 


The  values  of  tlie  parameters  B2  and  1)2  are  obtained  from  next 
matrix  equation  (§  3); 


K  ' 

1  ^9 »« 

Ba 

A 

^1. 

'b2i 

Ct 

^3» 

)fkl 

1>M 


The  necessity  of  converqency  needs  aqain  some  more  computation. 


A 

Asiy^(^-Ma](4-vi0 

c'lxAai  - 

C  .a 

T)”! 

A 

^  \)A»-  Ba 

%  Xi 

\  c\t- 

C,2 

L 

-  Cj, 


V  V  V  J 

-V-*»  ^  «-! 

J 

s 

.  ? 

■  .  ^  "  W  "*  - 

*  “  it*'  ■  ■  fc  •  « 

• '  sr%  r-*- 

•  •*/ .■  ■* 

'  •»'  .  .  * 

-  v'v  !-  >; 

.. 

'/•'a 

•"T.  •  -  •  •  ' 

vM'-a-:- 

! 

'.V  V.'.\  • 

V-‘'  -'A*A’ 

'a*  '•*  /**-*•* 
n‘a* 


^  b'„  V  cL''-'Vu.k«J  !»■ 

_  , ,  c“i'<-’Vw  b'., 

,2t-y-«)-j,l-/-V^  1  r  V'.i.b,,  -  bV  biij 

—  e.  .«,  o*  7\  L  ' 

4  j_l,'„.W.-  b'j.  1»0 

4  .-■'-''■’'''o:,*  b'*>  i“] 

=  _  -^' fv-x)  ^.. 

C-i:.  ^  ,  , 

The  positive  exponent,  e^^  included  in  and  has  disappeared. 
c.'^a  = 

V',2.U,  -  b'..  W] 

4  ’[V„,U,  -  W2,.v«] 

,  _-(v-’''' I  b«.b4,  - 


-v’C-'K-, 
•/  •.*  */  %  ^ 
:>  vA' 


r^-v --N  •^•v; 


.  ►  A  •  •  *  • . 


4  .  \  %  «‘#  -*J 


y'3.?3 


■'ai.  - 


_-,,b,-y-)  J.tY-')  8,a* 

.  e-W-''^ 


,  )  -(.'<'’‘)-ri(y-x) « 

cx\,.hx^  +  ^ 


The  positive  exponent,  e^2^^  included  in  and  b^^*  disappeared. 

The  numerators  of  B2  and  D2  contain  only  neoative  exponents  because  of 
the  presence  of  the  factors  a^^,  322*  3^3  ^24' 

4.3.  Values  of  the  parameters  A^,  Bp  Cp  Dp 


The  values  of  the  parameters  and  are  immediately  obtained  from 
the  matrix  equation  in  §  2.3. 


4Rn  Bz  ^ 


A,» 
C,  = 


The  values  of  the  parameters  and  are  obtained  from  the  boundary 
conditions  at  the  surface 

B,»  _ f 

S,~  ^  *- 

^  ^  \  ^  ^ A.r’'-  ^  1 

"  j-,-  ^  i- 


5.  Determination  of  the  stresses  and  the  displacements 


5.1.  Mathematical  procedure 


The  stresses  and  the  displacements  are  deduced  from  following  relations 
wherein  the  notations  of  §  1.  are  utilized: 


f 


-») 

■»• 

-ty>(z-H,.,  ) 

Di  e 

^  -b,  (‘-’'■-'ij 

oUaa. 

[A: 

Hi.,  ) 

4- 

>t_  -vwj.-  6t- 

Hi-.  )1 

v'l-W 

J 

tnr  ^ 

-Lra.-V  r  H 

I  E;  wi  L  .  ..  . 


J  yvx  L  /\i> 

^  .  a.  N 


-vnJ; 


The  stresses  in  cartesian  coordinates  are  calculated  as  explained  in 
appendix  1. 


5.2.  Stresses  and  displacements  at  the  surface  (z  =  0). 


=  O 


^  ^  lo 

^  Cv»i->^0C^»-')  *'* 


_^o^.  _ 1  TJvv.rV3s(v^«^')  [A,e^  s.r,e.^'’'3c*-^ 


!"*  ;y ,  (w^(*.■^ .  [To 
^  (^v*,vV0C^’-’')  Jo 

_  b  o. .  — It,  i.,.^')!  5;i~r^  -  iife^3l'''^’*''  ‘'' 


Y:^  X,  c**^'')  ■  C>*^  = 


I  YA 

\  Va 


(r;^<^) 


rj,f»^r).^.(v-‘>)  -.  ( VA 


^•C^Vv'.> 


'-V»‘I\.‘-jk' 

h\V*.‘^-'lCji 

C*. 

•-’>y-  .‘-I-"' 
‘-“v-yv' 

-/A  w‘* 

-/  •/  A 


■  v' 

■.•-■.'■\  •'.'t 

•..  %  N  % 

.*  •* 


5.3.  Stresses  and  displacements  in  the  first  layer  (0<h<Hp 

I 

[/I  s\^  ^  -r,vvnCW,-V'')\ 

+  1  (wirVTxG^^^)  +  C,e  J 

(vw^rVT.Cvw'^'l  I 

+  o, -e-  _  iJ-.Cie,  *-  J 

rr»rs  -  -  '  T.C'^'‘l  T',Cv«'«'l 

I  ±kl.C,  Z^''^  f ±^1-1^,  e  CA'- 

v^,-y,  v^.-yA,  J 

-  .  V*-  rT.M.T.C'""!  I  5,C~''-  «'■"'*''] ' 


■e 


■'  L  ».-)-• 

V),-y-i 

nth 

!*«  .\  V-  L  \*'  ((>.^-x')^^ 

n,  ^  '■  '  -n-  /JV.V*%r’-2^r)(V,VcASr’-2^r) 


-.^>*’1 1 


7«[|V)r]  .T,(W>^y  A  '  oUiA.  a 


.  it  f I  L^.lw^  - .  c,  -'■'3 


-i^tW,+-^.r<') 


*-  vii'-y** 

;h 

<0  '’^  \  -\ 

^  \  T,L‘^^)T«L‘Mr)  ^-w»V  >»Vi  J-^\y\ - 7,  \ 

[ — i;;? - 

tj'-  r  pt*v-f'--2Ar6»A<^ 

A  t  .  _Vm(^W,-V>)  -VV>V  rs  -W>ri  (Hi-I')  -  T) 


=  «  :n-  T.Cw>0 

Cj‘i-’1  .. 


- 


il  r  ^  ->-(H.-M  j.  r  ^->"'■l'^■■‘''^  cli~ 

V 


-  l^L  1  '*?,  Iv^r^  .T,  (>-M  [  >■)  A  ,e  “•  '"''''2  1 1,  A, f  '^  ■'  ''"'^ 


^  Q 

wrT,c-r).j,fw.^)c'’'“''cA^=  kri  r 


x_ 


.  -  i-  ,  •  -  .%  .*%-  .'  4.^  .V  _'^  ."V  .’^ 


with 

\  T,6*^r). 
K  Wi 


3,tw»f).  Ji(vhk) 
W\ 


) 

to’ 


5.4.  Stresses  and  displacements  in  the  other  layers. 

The  stresses  and  displacements  are  immediately  deduced  from  the 
relations  in  §  5.1. 

The  exponents  associated  to  the  parameters  are: 

-  in  the  second  layer  (H2^h/Hj+H2) 

-  in  the  third  layer  (Hj+H2  <'h<^Hj+H2+H2) 

^  ^vv,(K- H,- Hv) 

^  H24  H,- 

-  in  the  fourth  layer  (Hj+H2+H3  <^h<^Hj+H2+tl3+H^) 

V  -vn  (Vi- H,- Ha} 


APPENDIX  4 


Algebraical  Analysis  of  a  four-layer  anisotropic  System 
with  fixed  bottom,  full  slip  condition  at  the  first  and 
second  interface,  full  friction  at  the  third  interface. 


APPENDIX  4 


Algebraical  analysis  of  an  anisotropic  four  layered  structure  with 
fixed  bottom  and  full  slip  condition  at  the  first  and  second  interfaces. 


1.  Boundary  conditions. 


To  reduce  the  number  of  exponentials  we  write 
Ainry  ni 

i  r  ^  A  ■G- 

biM  r\,- ^  s  »• 

rs  -i  /■  ..  N  w.x;  -  •  Hi) 

Ci»vi  n;  riC'ni+bO®-  s'-' 

TV  .  .  .N--wtr.rHi+H,4-.-  Hi-4)  TV 


Further 

we  also  write 

V  - 

El 

Ey 

S'rfy'a- 

-i) 

“t  * 

y.  (jn^ 

E» 

W  = 

Ex 

Et, 

(*t3 
Jr  (rti- 

-O 

- 

L:. 

r 

E3 

Hj, 

?<  - 

rv.H. 

y 

wy  W  1 

MO 

/ 

•Z.  - 

.  w.  C >’• 

^  Mr  + 

Ma) 

)-■»  Wt  (  H 1 

+  Hx 

^3  + 

H.i,) 

where  H, ,  H,,  H-  and  H.  are  the  thicknesses  of  the  four  layers 


Boundary  conditions  at  the  surface  (z  =  0): 

A..*:’' >  H 

Xjis-o:  A.C''  -  b.  “ 

Boundary  conditions  at  the  first  interface  (z  =  Hj^); 

r»:  A,-*- c,  V  li, bx  +  Cl £■*’^^'^^4-  Pi 

Tix-  A,- b.e’' +  r.C, -^.b. 

W  *.  (-1+ ^')  At  ~  -  A  (lii+ es 

Boundary  conditions  at  the  second  interface  (z  =  Hj^  +  H2): 
fx  ••  Ai.  4-  Cx  b.  4-  p, 

t rz ..  A ,  -  bi  c-  V  T ,  Cx  -  rx  Pi  o 

VV!  [>4-*-V«0  At  -  ’’I  h  Cx  -  Si  Pit.  ^  a 

Boundary  conditions  at  the  third  interface  (z  »  +  H2  +  H^); 

W!  Cj -y  ^})  P)t  e 

L,  |(-^♦).A)Ai^e  +)«-«)  P*J 

U-.  (^k»>4V^s)A^  4-(»l»*h0  ^3^^*  4- ^'•  +  ^'3)^3  = 

Boundary  condition  at  the  bottom  (z  »  Hj  +  H2  +  +  H^): 

h/i  (-I+Jaa)A>4  -f'^t-^ii) 

If  S-\  1  Cxt^o 

•*'***t!L‘^ 


If  S.  <  1  A.4  * 


2.  Expression  of  the  boundary  conditions  in  matrixforin. 


2.1.  At  the  third  interface. 

In  the  equations  at  the  third  interface,  A^,  or  C^,  is  replaced  by  its 
value  obtained  from  the  fixed  bottom  condition. 

We  write  the  conditions  at  the  third  interface  in  matrixform 


We  invert  M, 


Hs-  ChA 


.^'1  A 

l5  =  — 


,  ^  ^  ,  f»-y) 


Jj(x-y) 

-  )<•  -  *' 

If  S4>  1 


\  4-e, 


ni= 


-1  +  a. 


L,y*i.»iA4VA)L'^ J 

UC-4V**)  -1 


We  write  A^,  and  in  function  of  and  D^. 

A  £Ms6^0-^tO'^)  - ^ 

f^y  e - - 

%3^CA~y\i) 

^  2r\^l20.v\(,C^^).^A  A-  ^  j^-y) 

)Dy  —  ’  *”' 

=  - 

^  :SLM<r(^0.  ^ (‘O 

Pb  =  - - ®- 

where  Mg(i,j)  are  the  constants  in 
We  write 

Tj^  - 

=  SM6(iOH462> 

SO  that 


IT" 

T"  .  ' 

A 

^31 

T« 

h 

u 


A4.5 


2.2.  At  the  surface 

Adding  and  substracting  the  surface  conditions  we  obtain 


lb,  = 

2.3.  At  the  first  interface. 


We  add  and  substract  the  first  two  conditions 

a  I.A  ij 


We  replace  Aj  and  Bj  by  their  values  obtained  from  the  surface  conditions 

(-itt.'jC,  ’'"j  ♦  t.i-r.)3>.[e:’'-  =  MO-'-’' 


We  solve  the  system 


^  fv  . 

v|_ls,c’‘ -  '  '  J  )■ 

The  positive  exponent  e^  has  disappeared. 


V. 


V^=  ‘iy.re"’’'-  (^+ 


,r  -3x 

‘  t-'- 'O’ [-'»«- ■'■  J 


For  wsao 


We  transform  the  w-cogdition  utilizing  the  t  ^^-conditions 

C,  -  r,  T,  Cic.  —  >k(Wi-0 

F  1  F» 

y.Cvi.-O 

=  FC^fT'-fv-'L  -Ft), 

We  replace  Cj  and  Dj  by  their  values 


For  mos 


9-a  =  >  ^z.= 


|[Aa].8,  TCie 


’  1  y. 

[Al]  -  11  D.]-- 

Writing  R'j®  ,  we  obtain  the  system 

R-1 

v^  -  1?3T>2.= 

n-i 

-  V  4  J.C,  -nP.  =  . 

and  by  adding  and  substracting 


o  ^  I 


2.4.  At  the  second  interface. 

We  add  and  substract  the  first  two  conditions 

=  [A3] 

a  (-i-rx)Ci  <-  ^  [A  3] 

We  replace  A2  and  B2  by  their  values  obtained  from  the  first  interface 
conditions  ..  ^ 

C.[C— 

We  solve  the  system 

'  .  .  \  ^  I  .  -.\n  t  A 


(Y-") 


t  e 


^  IXj  -  «-  Jl - 

‘  lY-«' 


Vz 


The  positive  exponent  has  disappeared. 

^or  yy>sM  ,  Vt*  ^  ^'v') 


►  *'•  ji 

■•■>■.•".■-■.••■5 


'•■'3 
-■-  -'•  .'*3 


ap' ' 


-v  •.  •.  •■  ' 

- ■  •  -.‘'l 

.1 


We  transform  the  v/-condition  utilizing  the t  ^^-conditions 

=  Vc,  ^3]] 

W=  k. 

We  replace  C2  and  D^  by  their  values 


Cx- 


iL-AJlLC'-^r.-KO 


-1J2('J-V) 

fc.  + 


(^-  Tj  ^  H 


X  ac’‘  ''•"i  c''’"'’-  2 

R.L 


-IN-’')  _2J^('^-’')-\ 
e  t  J 

2Ji  t.  «-  J  j  • 


T<»r  yvT«  aJ  ,  <S)\=  (>i- Jr-*- )1>>  > 

Writing  iPxs  ,  we  obtain  the  system 

fil, 


^^_\a-v)  _  ^  «  r»»T)3 

and  by  adding  and  substracting 

4  (-A4  j-i-ifci)  ('I- 5>  *-03)5)3 

“ibj  4  ( -A-S;,- fil3)C-3  4 


(S>, 

-  ik 
li. 


[>A4  ^3  4  ®  i)  T)>=. 


3.  Resolution  of  the  system  of  boundary  conditions. 

We  have  from  the  boundary  conditions  at  the  third  interface  (  §  2.1) 

— -  T  pM-ft*  +  Pa- ^■‘3 

_ »  Pw 0*3 1- 

a  ^  r  Pi.-b* 

_ IH, .  b*  +  p*-  P‘3  - 

lt3(A-V»3l 

We  replace  A^,  and  in  the  last  equations  of  §  2.4 


We  solve  the  system  in  and 


M.IO 


£)»  L 

■n  -(Tt-vl  .  o  -fi-y)! 

4-  S-b  '•‘42*''^  *  *-  -  ■*>  )•  ”1 

V3 

The  positive  exponents  e^^  e^3^^  and  e^^  have 

disappeared. 

T)a  =  2,,C''-‘'^0  [W.  f..  e'-'l"-''’ 

. . ^ 

^  (yn  4  P„  Pai  -  PnP*0 

/  «  o  \ 

4-  (P..  -  P'"-  ‘‘>0 

w  «,= «  v„  =  L(^- ’*■>  *  -  ’’’■  ^'‘’) 

The  numerators  of  B^  and  D^  tend  both  to  zero  and  the  denominator 
tends  to  a  constant  value. 


4.  Values  of  the  parameters  . 

4.1.  Values  of  the  parameters  A^,  C^,  0^. 

The  values  of  the  parameters  A^,  B^.  C^.and  are  obtained  from  the 
boundary  conditions  at  the  third  interface  in  which  B^  and 
are  replaced  by  their  values  from  §  3. 

s  ^  Pii-Pai  - 

^  Sylhvh^-  Pn-PAi) 


2)3.  V  Pa.  64  +  Pii 

^  7Si(4-y^i) 

.t.  r  ft  O  \  -nl-x-y) 

=  ^  I  (Pia-Ta.  -  e-  •  «- 

_  Pii  ''^ 

■  6,  .  .T  1 


^.7n(»-v) .  ,  Pi,.  Pia)]  V, 


The  positive  exponent  e'^'P)  has  disappeared.  Although  the  presence 
of  the  constant  S2(p22^41  “  ^21^42^  numerator  converges  to  zero 
because  of  the  factor  Q2. 


C35  ^ 0 


,  I  »a,.p„  -  Pal.  P.o  [P«  Pi,- 


m 


=  ^  y  L^a>  Pa  -  P.>-  '^^^)  V  (piil-  -  ^Av 

L 

4.2.  Values  of  the  parameters  A2,  B2,  C2,  D2. 

The  values  of  C2  and  D2  are  obtained  from  the  relations  established  in  §  2.4. 


C, 


Vr 


T).=  ^ 


V-! 


R.)i’lv-')  .  [.-H*  R.)J. 

V  6i  •■  CTj  i'>l'‘''>  +  P»3 


f  Ei 
8 


The  numerator  in  D2  converges  because  of  the  presence  of  the  factor  R2. 
The  value  of  obtained  from  the  relation  established  in  §  2.4. 

The  value  of  B2  is  obtained  from  the  last  relation  of  §  2.3. 


4.3.  Values  of  the  parameters  A 


The  values  of  the  parameters  and  Cj  are  obtained  from  the  relations 
established  in  §  2.3. 

^  4  I>»] 

* 

For  the  determination  of  Aj,  we  need  the  value  of  D^e”^!^. 


-fiX 


t-  f *  (4-j, )£'■”] 


The  values  of  and  0^  are  obtained  from  the  surface  conditions. 

B,  *  — - -  [  5,  -  J-O  -A.e- -  2J',C‘, 

5.--1  ^ 

3)  2  _  _ _  f  -A  -  -ly\\S>)  C,  t 

5.  Relations  for  the  stresses  and  the  displacements. 


The  relations  for  the  stresses  and  the  displacements  are  completely 
the  same  as  those  developped  in  appendix  3,  by  replacing  the  parameters 
A.,  [).  by  their  adequate  values. 

The  relation  for  the  verticale  displacement  is  again  undeterminated 
at  the  origin  (  m  =  0),  The  problem  is  solved  in  exactly  the  same 
way  as  developped  in  appendix  2  (§  4.). 


EXPLANATORY  NOTICE 


I 


ISIC  FOUFC- LAYERED  (FLxy)  USAE  Fage  1 

F‘our~-!l  aiyered  System  F’rogram  : 


This  pregram  is  civailable  in  two  versions: 

-•  EXECLIiADLE  version,  in  which  the  executable  progr  am 
is  made  up  of  only  one  block  with  automatic  loading. 
"  SOURCE  version,  in  which  all  controls,  data,  and 
all  modules  are  in  separateid  files. 

F 1  op p  y  cJ  i  s |:;  c  on  t  en  t  s  : 


In  following  text,  ;;  means  A  =  ANISOTROPIC 

I  =  ISOTROPIC 
y  means  P  =  PARTIAL 
S  =  SLIP 

Each  floppy  disk  contains  fallowing  main  files: 

AUTOEXEC.BAT  - >  automatic  program  lociding 

LOOO.BAT  - >  display  of  introduction  logo 

FLxyLO.TXT  - >  introduction  text  on  screen 

FLxyNO.LOG  - >  introduction  text  for  printer 

FLx yMO . TXT  - >  this  not ice 

$ EXECUTABLE  floppy  disk 

This  disk  contains  the  fallowing  files  in  addition  to 
the  main  filcjs: 

FLx y.  EXE"  - >  executable  pr"ogram 

FLxy.DAT  - >  data  file  for’  demonstration 

FLxy.LST  - >  result  file  for  demonstration 

* SOURCE  floppy  disk 

This  floppy  disk  contains  the  following  files  in  additicjn  in 
the  main  files: 

FI..,;:  y VE I  i  . .  >  revi  si  on 

F L.. y  1 . 1  OR  - >  m  i  n  m o d  u  1  e 

FLx  y 2 . 1  T.1R  - >  subr out  i  nes  DOMEC 

DOERA 

FOCAL 

PC)  I  NT 

PAS 

CHECH 

ERROR 

FLxy3.F0R'  - >  subroutines  ECHDE 

V I N I T 

ZERO 

FINIT 

F  L.  y  '1 .  f”  C!)  R  - >  s  i.i  b  r  o  u  t  i  n  t?  s  P  4  ^  4  2 

PCT22 
PCT42 
S0N42 
SDN22 
P4444 
P244  2 
CONST 


»•.  r. 

.  •'Jl  •  •  •  •  ' 
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)|( ))( It; «  i  « Ik  i|( « )tc  )i(  i|t  )(c « Ik  « )ti  It:  )|c  Ik !(:  >k  !|c  Ik  )t: « Ik  i|(  Ik  «  « t  !|(  i|(  4:  :|c  It!  i|(  i|c  i|( ))!  It:  i|(  lit i|<  lit  Ik  i|( « Ik  !|t  *  i|(  !|(  i|(  i|(  )|(  ]|(  iti  )|(  i|( 


FLityS.FQR 


FLity^.FOR 

FLityy.FQR 


FLitya.FQR 

FLi:y9,F0R 


Con  f i gur  at i an i 


—  >  Bubr  OLiti  nes 


— >  subroutines 


iubrouti nes 


subroutines 

subroutines 


BES  J  1 

BJ0J2 

SURFA 

OOIJC 1 

CrJUC2 

CDUC3 

C0LIC4 

T I TRE 

MOD  IF 

BIMF^S 

COMSU 

FONG 

CCDLIl 

FONG  1 

F0NG2 

FDNG3 

F0NG4 

FONGS 

F0NG6 

EGHEF 

IMDON 

IMRES 

MENU 

EGDON 

LEDON 

AFDMG 

AFDTR 

AFPOG 


A  v'y  v>.t!3 


m. 


lEtM-PG  (GjXTjAT)  with  at  least  2561<:Et,  1  ou  2  diskette  drives, 

00  col.  screen  (monochrome  or  color),  math,  coprocessor, 
matr  i ;;  pr  i  nter  . 


Running  of  the  program; 


In  the  next  presentation  <ENTER>  means  action  of  key  <- 
Insert  the  EXEGUTABLE  disk  into  drive  B. 

Locid  DOS,  if  necessary,  then  type  DIR  Ei;  and  PATH  A:  \ 
FLB>  FLxy  < ENTER) 


Printing  of  the  results; 


on  screen 
on  printer 


:  FLB>  type  FLxy.LST  < ENTER) 

;  FLB>  type  FLxy.LST  >prn  <ENTEF<> 


Gonsul  ti^ti on  of  this  notice; 


on  screen 
on  printer 


FLB>  type  FLxyNO.TXT  < ENTER) 
FLB)  type  FLxyNO.TXT  )prn  CENTER: 


•-  .•> 

•  •  >.*  O  V 
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Disk,  preparation  (or  normal  running: 


Format  ai  system  disk  with  COMHAND.COM,  and  files  needed  for 
r Linn i  ng  AUTOEXEC .  EiAT . 

Preparation  of  CCDNFIG.SYS  with  files=10,  device=ansi,  buffers-lf< 


Alterations  in  source  files: 


By  teiit  editor  EDLIN  or  any  other  text  editor,  program 
statements,  may  be  altered.  ‘  * 

Commands  for  EDLIN  are: 

FLB  >  ed  1  i  n  x  >;  x  x  n .  f  or  < ENTER  > 

where - >  file  n ame  to  be  mod i f i ed  < x  x x x  n . ( or ) 


nD  — > 

n  —  > 

nl,n2L  — > 
E  '  — > 


erase  line  number  n 

displays  line  n  for  alteration,  type  the 
correct  statement. 

displays  lines  between  number  nl  and  number  n2 
ends  session  and  retuirns  to  DOS 


Compi  1  £iti  on : 


After  alteration,  the  new  version  of  the  module  has  to  be 
compi 1 ed . 

Insert  Professional  fortran  compiler  into  drive  A. 

FL.,A>  b:  <  ENTER  > 

FLB>  path  a:\  < ENTER > 

FLB>  profort  xxxx/li  < ENTER > 

Linkage: 


After  correct  compilation  has  taken  place,  a  new  executab  1 1:.; 
program  has  to  be  created. 

Insert  diskette  with  FORTRAN  libraries  into  drive  A. 

FLB  >  1  i  n  FLx  y  1  +FLx  y 2+ .  .  .  ,  FLx  y ,  CON :  ;  <  ENTER  > 

Note  : 


To  obtain  introduction  logo  on  printer 
FLB>  type  FLxyND.LOG  >PRN  <ENTER> 


u  ^1' 
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t$$*tt*tM***t*t**t$******1(********************************************** 

Presentation  of  different  possible  screens: 


'  -  '=. 


■■ ‘‘v 

.  .  •.  5  .  < 


FD-USAE-vers.  2.00  1986 

MAIN  MENU  scren  : 


Strains  and  stresses 


5^'. 


i  n  USAE-2 

a  four-layered  system. 

tt********************************************************************** 


1.  Data  retrieval  in  a  file 

2.  Data  saving  in  a  file 

3.  Screen  displaying  and/or  alteration  of  system  data 

4.  Screen  displaying  and/or  alteration  of  traffic  data 

5.  Screen  displaying  and/or  alteration  of  computation  coordinates 

6.  Input  of  intermediate  depths 

7.  Program  start 


Your  choice 


0  for  stop  —  : 


Screen  of  CHOICE  1; 


Strains  and  stresses 

i  n  USAE-2 

a  four-layered  system. 

*$*$$$**$**$*$**$$*$****$$$*$***$**$$$$*$$*$*$*******$*$$$*$***** 

Data  retrieval  in  a  file. 


A  FLxy.DAT  file  contains  base  data. 

The  user  can  define  another  file  whose  name  has  to  be  written  in 
B  characters  (format  XXXX.DAT). 


Name  of  chosen  file  or  FLxy.DAT 
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*********  **$*  ******  *******-^  *  ^^^^*^^^^* 

Screen  of  CHOICE  2: 


Strains  and  s  t  r  e  s  s  e  s 

in  IJSnL  2 

a  four-layered  system. 

************************************************************************* 
D  a  t  £<  s  a  V  i  n  g  in  a  f  i  1  e . 


A  f”L:;y.DAT  file  contains  base  data. 

The  user  can  define  another  file  whose  name  has  to  be  written  in 
8  charactcjrs  (format  XXXX.DAT). 

Name  of  chosen  file  or  FLxy.DAT  : 


Screen  of  CHOICE  3 


Strains  and  stresBt.‘s 


in  Ui:.)fU:'-2 

a  four-layered  system. 

************************************************************************* 


System  data. 


Layer 

Modul us 

Poisson’s  r. 

Thi  ckni 

1. 

400000 . 0 

0.  16 

20 . 0 

1?  • 

1 OOOOO . O 

0.25 

20.  0 

3 . 

1 0000 . 0 

0.  50 

30 . 0 

4. 

1 000 . 0 

0.  50 

9000. 0 

Friction  ra. 


1 . 0000 


1 . 000f» 
1 .  Of) 00 


For  return  =1 


or  a  1 1  er 
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*$***$********$********************************************************** 

Screen  o-f  CHOICE  As 


Strains  and  stresses 


USAE-2 


a  four-layered  system. 

*t$*$$*t*$*********$***t**tt$********t****t****************t********t*** 

Number  of  circular  loads  =  2 


Radius 


Pressure 


Dist.  H 


Dist.  y 


1 1 . 450 
1 1 . 450 


7.900 

7.900 


0.000 

34.350 


O.  000 
0 . 000 


For  continuation  =  1  or  alter.  =  2  - 


Screen  of  CHOICE  5: 


Strains  and  stresses 

in  USAE-2 

a  four-layered  system. 

t***t****$*t*t***$****t*******t*$t***$*$**t****$*$$$ti.**tt**************** 

Number  of  computation  coordinates  =  2 


0 .  CiOO 
17. 175 


O,  000 
0.  000 


For  continuation®  1  or  alter.  =  2  -  : 


Screen  of  CHOICE  63 


Strains  and  stresses 

in  USAE-2 

a  four-layered  system. 


Positions  of  stress  computations  in  depth  out  of  interfaces. 

(max  22) 
number  of  positions  i 


tor 
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tt%M  ******  **t***ttt**tt*t**t*****t****t$itt*t*tt*  it*****  t****t***.*t***  ******* 

Screen  af  CHOICE  7: 


Strains  and  stresses 

in  USAE-2 

a  tour-layered  system. 

*******************************************************************************.. 

initial  computation  interval  (0.1  is  ^eneraly  small  enough)  s 


Next  screen: 


Strains  and  stresses 

in  USAE-2 

a  tour-layered  system. 

********************************************* ************t*********************i 


choice  ot  scale  i 
***************** 

allowed  choices  : 

1  -  thickness  of  Irst  layer 

2  -  thickness  of  2  first  layers 

3  -  thickness  of  3  first  layers 

4  -  thickness  of  4  layers 

5  -  load  radius 


suggested  solution  :  S 

your  choice  : 


Screen  for  execution: 


Strains  and  stresses 

in  USAE-2 

a  four-layered  system. 

******************************************************************************* 


computation  start .  be  patient!! 


be  even  more  patient...!! 


•  ‘  »  • 

•  v'  •  r  • .  . 

"  •>  'j-  ' 

\  O  . 


n 

.  *  e  ^  t  ' 

,  «\V. 


»  A  »• 
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'*«»»»»**«»»»« 44*»*»»* «»***«»»**»*»»«***»«» 

Screen  For  results; 


Strains  and  stresses 

in  USAE-L: 

■  a  Four -1  ayereiJ  ‘..ysleiii. 

:t.**i.*t.t$**t  $»**.*t**$*ii.t*****t**-****tt****ttt****t*********  ***■*■***»■****  ***  »*■*■**** 

I 

Ia  File  FL>:y.LST  contains  base  results. 

I  The  user  can  deFine  another  File  Mhose  name  has  to  be  written  in 
'b  characters  (Format  XXXX.DAT). 

! 

■  Name  oF  chosen  File  or  FLi;y.LST  : 
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'ample  of  results  (FLxy.LST) 


our-layerea  system  program  isotropic  -  partial  triction 

U5h£-2  dtpirtikcnt  iis  constructius 

isic 

av.dt  rhopital,  27  h 
7060  aoni  Odtiout 

iiiuiuiiuiiiiiiiuitttiiniiiiiiuiiiiiiiiitiiiiiuiiiitiiiiitiniiiniitiiiiiiiiiiiiitiiniiiuiiiiiiiiiiiiui 

lechanical  data 

«  4  » I »  «  4  a  »  *  a « y  0 


Yuiiiiij'u  iiiuilulus 

F.  '  k  r  at  1 u 

t  111  Cl.llUbS 

1  r  1  c  1 1  nil  1 

400000.0 

0.16 

20.000 

1 . 00 

100000. 0 

0.25 

20.000 

1 . 00 

10000.0 

0. 50 

30.000 

1 . 00 

1000. 0 

O.  50 

9000.000 

ra-ffic  data 

aa********* 

load  radius  pressure 

y 

1  11.450  7.900 

0.  000 

0 . 000 

11.450  7 . 900  34 . 350  0 . 000 


ISIC  FOUR-LAYERED  (FLxy)  USAE  Rage  10 

tittutitittiuititttttuttititttttttttittutttlittttittlitttittulitttitit 

position  1  K  =  0.000  y  =  0.000 


depth 

layer 

SI 

»y 

sz 

tyz 

tiz 

txy 

0.000 

1  tt 

10.I6&0 

■11.&793 

7.9000 

0.0000 

20.000 

1  It 

-3.4792 

-4.5018 

1.4988 

0.0000 

-0.4167 

0.0000 

20.000 

2  11 

-0.5790 

-0.8163 

1.4989 

0.0000 

-0.4166 

0.0000 

40.000 

2  It 

-2.2733 

-2.5873 

0.1873 

0.0000 

-0.0611 

0.0000 

40.000 

3  11 

-0.173i 

-0.1997 

0.1873 

0.0000 

-0.0611 

0.0000 

70.000 

3  It 

-0.3907 

-0.4035 

0.0540 

0.0000 

70.000 

4  It 

0.009S 

0.0082 

0.0540 

0.0000 

-0.0046 

0.0000 

9070.000 

4  It 

0.0000 

6.0000 

0.0000 

0.0000 

depth 

layer 

si 

s2 

s3 

epsl 

eps2 

eps3 

0.000 

1  It 

11.6793 

10.1660 

7.9000 

0.2197E-04 

0.1758E-04 

O.llOlE-04 

20.000 

1  tt 

1.5335 

-3.5138 

-4.5018 

0.7040E-05 

-.7597E-05 

-.1046E-04 

20.000 

2  tt 

1.5793 

-0.6595 

-0.8163 

0.1948E-04 

-.8502E-05 

-.1046E-04 

40.000 

2  tt 

0.1888 

-2.2749 

-2.5873 

0.I404E-04 

-.I675E-04 

-.2066E-04 

40.000 

3  It 

0.1973 

-0.1836 

-0.1997 

0.3890E-04 

-. 1824E-04 

-.2066E-04 

70.000 

3  It 

0.0540 

-0.3908 

-0.4035 

0.4S12E-04 

-.2I60E-04 

-.2352E-04 

70.000 

4  It 

0.0545 

0.0091 

0.00S2 

0.4582E-04 

-.2230E-04 

-.2352E-04 

9070,000 

4  It 

0.0000 

0.0000 

0.0000 

O.OOOOE^OO 

O.OOOOE^OO 

O.OOOOE^OO 

depth  1 

layer  ulil 

v(yl 

m(z) 

ei 

ey 

ez 

0.000 

1  tt  -.2184E-03 

-.1909E-10 

-.9046E-02 

0.1758E-04 

0.2197E-04 

O.llOlE-04 

20.000 

1  It  0.988SE-04 

0.8642E-11 

-.8835E-02 

-.7497E-05 

-.1046E-04 

0.6939E-05 

20.000 

2  It  0.988SE-04 

0.8642E-11 

-.8835E-02 

-.7497E-05 

-.1046E-04 

0. I848E-04 

40.000 

2  tt  0.3017E-03 

0.2637E-10 

-.8561E-02 

-.1673E-04 

-.2066E-04 

0.1402E-04 

40,000 

3  tt  0.3017E-03 

0.2637E-10 

-.8561E-02 

-.1673E-04 

-.2066E-04 

0.3739E-04 

70.000 

3  tt  0.3844E-03 

0.3360E-10 

-.7432E-02 

-.2159E-04 

-.2352E-04 

0.4S11E-04 

70.000 

4  It  0.3844E-03 

0.3360E-10 

-.7432E-02 

-.2t59E-04 

-.2352E-04 

0.4511E-04 

9070.000 

4  It  O.OOOOE^OO 

O.OOOOEtOO 

O.OOOOEeOO 

O.OOOOEtOO 

O.OOOOEtOO 

O.OOOOEtOO 
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Page 
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tItIttllltItiltlltlItttMlllllltlltMllltIUtttItttttttlItttUttItttttlltll 

position 

2  X 

=  17.175 

Y  = 

0.000 

depth  I 

layer 

sx 

sy 

sz 

tyz 

0.000 

1  U 

3.0499 

B.82S9 

0.0000 

0.0000 

20.000 

1  It 

-2.0093 

-4.3575 

1.1255 

0.0000 

20.000 

2  It 

-0.3243 

-0.8505 

1.1257 

0.0000 

40.000 

2  It 

-2.4172 

-2.7155 

0.1997 

0.0000 

40.000 

3  It 

-0.1028 

-0.2077 

0. 1997 

0.0000 

70.000 

3  It 

-0.4084 

-0.41&& 

0.0554 

0.0000 

70.000 

4  11 

0.0090 

0.0082 

0.0554 

0.0000 

9070.000 

4  It 

0.0000 

0.0000 

0.0000 

0.0000 

0.000 

1  It  0.8259 

3.0499 

0.0000 

0.20B4E-04 

0.4094E-05 

-.47508-05 

20.000 

1  It  1.1255 

-2.0893 

-4.3575 

0.5392E-05 

-.3930E-05 

-.10518-04 

20.000 

2  It  1.1257 

-0.3243 

-0.8505 

0.1419E-04 

-.3931E-05 

-.1051E-04 

40.000 

2  It  0.1997 

-2.4172 

-2.7155 

0.14B3E-04 

-.178BE-04 

-.2161E-04 

40.000 

3  It  0.1997 

-0.1828 

-0.2077 

0.3949E-04 

-.17B8E-04 

-.2161E-04 

•  k.*V  m'*  b 

70.000 

3  It  0.0554 

-0.4084 

-0.4166 

0.4679E-04 

-.227BE-04 

-.24018-04 

70.000 

4  It  0.0554 

0.0090 

0.0082 

0.4679E-04 

-.2278E-04 

-.2401E-04 

9070.000 

4  It  0.0000 

0.0000 

0.0000 

0. 00008*00 

0.00008*00 

O.OOOOE*00 

f  1 

ipth  layer  u(x) 

v(yl 

w(z) 

ex 

ey 

ez 

0.000 

1  It  0.2272E-18 

-.1565E-10 

-.8979E-02 

0.4094E-05 

0.20848-04 

-.4750E-05 

20.000 

1  II  -.4354E-19 

0.7889E-11 

-.8961E-02 

-.3930E-05 

-.10518-04 

0.53928-05 

20.000 

2  It  -.3440E-1B 

0.78a9E-ll 

-.a96IE-02 

-.3931E-05 

-.10518-04 

0.1419E-04 

40.000 

2  It  -.7092E-18 

0.1622E-10 

-.86988-02 

-.178BE-04 

-.21618-04 

0.14838-04 

40.000 

3  II  -.7092E-18 

0.1A22E-10 

-.8698E-02 

-.17BBE-04 

-.21618-04 

0.39498-04 

70.000 

3  It  -.7879E-I8 

0.1803E-10 

-.7520E-02 

-.227BE-04 

-.24018-04 

0.46798-04 

■  1  .  .  ,  . 

70.000 

4  It  -,7879E-18 

0.1803E-10 

-.7520E-02 

-.227BE-04 

-.24018-04 

0.4679E-04 

■\  *  ■'  ■*  '*•  ] 

9070.000 

4  It  O.OOOOE^OO 

O.OOOOEtOO 

O.OOOOEiOO 

0.0000E*00 

0.00008*00 

0.00008*00 

'  -w'  v'  V  1 

'•'."■'.■'•'AW 

utilized  symbols: 


sx  :  normal  stress  in  the  x-direction 

sy  :  normal  stress  in  the  y-dl recti  on 

sz  ;  normal  stress  in  the  z-direction 

tyz  :  shear  stress  in  the  yz  plane,  parallel  to  y  or  z 

txz  :  shear  stress  in  the  xz  plane,  parallel  to  x  or  z 

txy  :  shear  stress  in  the  xy  plane,  parallel  to  x  or  y 

si  :  maximum  principal  stress 

s2  :  medium  principal  stress 

s3  :  minimum  principal  stress 

epsl  :  principal  strain 

eps2  :  principal  strain 

eps3  :  principal  strain 

u(x)  :  displacement  in  the  x-direction 

v(y)  :  displacement  in  the  y-direction 

w(z)  :  displacement  in  the  z-direction 

ex  :  strain  in  the  x-direction 

ey  :  strain  in  the  y-direction 

ez  :  strain  in  the  z-direction 

The  normal  stresses  are  taken  positive  when  they  produce  compress 

and  negative  when  they  produce  tension. 


